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PREFACE 


The  purpose  of  this  study  was  to  determine  if  it  were  possible  to  construct  a  theoretically 
based  model  that  could  produce  ''reasonable"  vertical  profiles  of  drop  size  distributions 
and  associated  atmospheric  parameters  in  very  low  stratus  clouds  and  their  associated 
subcloud  regions  and  at  the  same  time  have  a  complexity  somewhere  between  the 
“simple”  empirical  model  currently  used  by  the  U.S.  Army  and  highly  complex  academic 
level  models.  Moreover,  the  inputs  required  to  solve  the  highly  nonlinear  equations  that 
characterize  the  model  would  consist  of  a  small  number  of  ground  level  (2  m)  conven- 
tional  atmospheric  observations,  that  is,  no  vertical  profile  measurements  would  be  nec¬ 
essary.  Finally,  the  model  design  should  be  simple  enough  to  be  executable  on  a 
comparatively  small  computer. 

The  overall  approach  (in  an  attempt  to  accomplish  the  above)  is  nontraditional  for  very 
low  stratus  clouds  and,  consequently,  has  become  controversial.  There  are  those  who 
feel  that  the  physics  and  thermodynamics  of  rising  clusters  of  droplets  are  not  appropriate 
for  the  vertical  structuring  of  low  stratus  clouds.  However,  results  to  date  from  the  model 
are  so  encouraging  that  I  feel  it  must  be  a  good  analog  for  constructing  the  vertical 
structure  of  very  low  stratus  clouds. 

There  are  others,  including  the  coauthor  of  this  report,  who  insist  that  I  am  trying  to  get 
too  much  out  of  the  model;  it  has  been  overextended.  These  extensions  are  driven  by 
application  needs,  and  I  feel  they  will  provide  the  ballpark  estimates  required. 

As  for  the  equations  used  in  the  model,  some  are  well  founded  in  the  literature,  others 
are  nothing  more  than  modern  updates  of  those  used  in  the  past,  and  some  are  of  our 
design.  Finally,  some  ideas  and  constraints  are  based  on  my  interpretations  of  relatively 
recent  approaches  and  results  by  other  authors. 

Henry  Rachele 
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SYMBOLS 


a  =  surface  tension  constant  =  0.153  dynes/cm  K  (in  equation  (15)) 

a  =  temperature  radiation  constant  (when  used  in  lapse  rate  equation) 

aj  =  constant  coefficient  in  polynomial 

A  =  a  double  exponential  constant  (in  equation  (2)) 

A  =  coefficient  defined  by  equation  (55a)  (except  in  equation  (2)) 

Ac  =  net  radiation  effect  in  the  cloud  =  an  unknown  function  of  z 

Aw  =  activity  of  water  (Aw  =  fr  when  8=1)  =  exp  (-*?/(Mw/Mcj)) 

Aw,  =  threshold  water  activity 

B  =  a  double  exponential  constant  (in  equation  (2)) 

B  =  coefficient  defined  by  equation  (55b)  (except  in  equation  (2)) 

b  =  surface  tension  constant  depending  on  water  soluble  substance 

(in  equation  (15)) 

b  =  an  adjustable  parameter  related  to  "vertical"  entrainment  from  layer 

above  cloud  top  (when  used  in  lapse  rate  equation) 
c  =  an  adjustable  parameter  related  to  "vertical"  entrainment  from  layer 

above  cloud  top  (when  used  in  lapse  rate  equation) 

C  =  a  double  exponential  constant  (in  equation  (2)) 

C  =  lumped  value  defined  by  equation  (55c)  (except  in  equation  (2)) 

Cp  =  specific  heat  of  dry  air  at  constant  pressure 

Cpc  =  specific  heat  of  water  vapor  at  constant  pressure 

C'p  =  weighted  sum  of  specific  heats  as  defined  by  equation  (22f) 

D  -  lumped  quantity  defined  by  equations  (27c)  and  (55d) 

Dv  =  diffusion  coefficient 

E  =  Contribution  of  "horizontal"  entrainment  to  dT/dz  in  cluster  sphere 

(in  equations  (22a)  and  (22b)) 

E  =  force  necessary  to  accelerate  the  mass  being  entrained 

per  unit  of  total  mass  in  sphere  (in  equations  (35)  and  (38)) 
es  =  partial  pressure  of  water  vapor  at  saturation 

ev  =  partial  pressure  of  water  vapor 

Fa  =  drag  force  per  gram  of  mass  in  sphere 
f  =  relative  humidity  of  the  ambient  air 

fr  =  relative  humidity  of  the  air  at  the  droplet  surface 

fr  =  relative  humidity  of  the  ambient  air  at  (z  -  Az) 

fr<  =  relative  humidity  of  the  air  at  the  droplet  surface  at  (z  -  Az) 

g  =  acceleration  due  to  gravity 

G  =  3(ln  fr)/3r  when  fr  is  defined  by  a  modified  form  of  a  Hanel  droplet  growth 

equation 

H  =  -  3(ln  fr)/3T  when  fr  is  defined  by  a  modified  form  of  a  Hanel  droplet  growth 

equation 

K  =  coefficient  of  heat  conduction  or  exchange  coefficient 

ki  =  mixing  rate  for  momentum 
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=  mixing  rate  for  heat 
=  specific  enthalpy  of  vaporization  of  water 
=  mass  of  sphere 
=  mass  of  dry  particle  material 
=  mass  of  solute 
=  molar  mass  of  solute 
=  mass  of  water 
=  molar  mass  of  water 

=  mass  of  water  condensed  on  the  particle  per  unit  mass  of  dry  material 
=  number  density  of  particles  of  radius  r 
=  the  cumulative  number  density  of  particles  of  radius  r 
=  concentration  of  cloud  condensation  nuclei 

=  number  of  droplets  of  size  class  i  per  cubic  centimeter  (i  =  1,2,..  .) 
(except  in  equation  (53)) 

=  the  number  density  corresponding  with  mode  i  of  a  bimodal  distribution 
(in  equation  (53)) 

=  Nj  (number  of  droplets  of  size  class  i  per  cubic  centimeter)  at  the  reference 
height 

=  number  of  moles  of  dry  solute 
=  number  of  droplet  size  classes 
=  mole  number  of  water 
=  Pd  +  ev  =  total  pressure 
=  partial  pressure  of  dry  air 
=  specific  humidity 
=  droplet  radius 

=  initial  dry  particle  radius,  r(f  =  0) 

=  r^  for  size  class  i 

=  radius  of  droplet  in  size  class  i  (i  =  1,2,...)  (except  in  equation  (53)) 

=  the  mode  radius  for  mode  i  of  a  bimodal  distribution  (in  equation  (53)) 

=  radius  of  cluster  sphere 
=  specific  gas  constant  for  water  vapor 
=  specific  gas  constant  for  dry  air 
=  specific  heat  of  liquid  water 
=  specific  heat  of  solid  or  dissolved  substance  i 
=  supersaturation 
=  time 

=  temperature  of  cluster  in  °K 
=  temperature  of  the  ambient  air  in  °K 
=  standard  temperature  (273.16  K) 

=  temperature  of  the  clusler  in  -K  at  the  reference  height 
=  temperature  of  the  ambient  air  in  °K  at  the  reference  height 
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=  number  of  moles  of  ions  formed  from  one  mole  of  solute 

=  specific  volume  of  pure  water 

=  effective  ascent  rate  of  the  ascending  cluster 

=  reference  level  value  of  w 

=  vertical  acceleration  of  cluster 

=  grams  of  condensed  matter  per  gram  of  dry  air 

=  grams  of  liquid  water  per  gram  of  dry  air 

=  grams  of  solid  or  dissolved  substance  i  per  gram  of  dry  air 

=  grams  of  water  vapor  per  gram  of  dry  air  =  mixing  ratio  for  water  vapor 

=  expression  defined  by  equations  (33)  and  (55e) 

=  height  above  ground  level 
=  reference  height  (2  m) 

=  height  step  size  in  numerical  integration 
=  entrainment  parameter 

=  ratio  of  the  molecular  weight  of  water  to  the  mean  molecular  weight  of  dry  air 
=  0.621972 

=  exponential  mass  increase  coefficient 
=  value  of  r/  at  infinite  dilution 
=  potential  temperature  for  damp  air 
=  total  mass  in  sphere  divided  by  volume  of  sphere 
=  density  of  dry  air 

=  density  of  dry  air  at  reference  height 
*  density  of  dry  particle 
=  density  of  pure  water 
=  density  of  ambient  damp  air 
=  surface  tension  of  a  droplet 
=  extinction  coefficient 

=  the  standard  deviation  for  mode  i  of  a  bimodal  distribution 
=  surface  tension  of  pure  water  at  standard  temperature 
=  practical  osmotic  coefficient 
=  curvature  correction 

=  5  (AW))  =  curvature  correction  when  water  activity  is  at  the  threshold  value 
=  parameters  in  Hanel  and  Lehmann's  approximation  equation  for  the 
exponential  mass  increase  coefficient 


A  prime  (')  is  used  to  indicate  an  environmental  quantity  (i.e.,  outside  of  sphere), 
(a'  and  b'  are  exceptions  to  this  statement.) 


H 


1.  INTRODUCTION 


The  Army  has  a  need  to  assess  the  effects  of  the  atmosphere  on  the  performance  of 
modern  sophisticated  weapons  and  surveillance  systems  that  employ  electro-optical  (EO) 
sensors.  These  atmospheric  effects  assessments  are  used  by  people  involved  in  re¬ 
search,  analyses,  weapons  design,  wargaming,  deployment,  and  tactics.  Of  particular 
interest  is  the  ability  of  a  sensor  to  "see"  a  target  along  slant  paths  through  the  lowest 
kilometer  of  the  atmosphere,  especially  when  moisture  in  the  form  of  damp  haze,  fog,  or 
low  clouds  is  present.  To  assess  the  slant  path  effects,  one  requires  knowledge  of  the 
vertical  variations  of  the  physical  quantities  of  these  moist  environments,  and  especially 
the  drop  size  distribution.  These  characteristics  are  currently  obtained  from  models  that 
are  based  on  field  data;  that  is,  they  are  empirical.  These  models  generally  satisfy  three 
user  preferences:  first  the  models  must  be  as  simple  as  possible,  but  produce 
"reasonable”  results;  second,  the  models  can  be  initialized  with  a  minimum  number  of 
surface  data  inputs,  and  third,  the  models  should  not  be  computer  intensive. 

The  purpose  of  this  study  was  to  develop  a  theoretically  based  model  to  augment  or  re¬ 
place  the  empirical  models  and  also  satisfy  the  user  needs  and  preferences  mentioned 
above.  In  particular,  we  wanted  a  model  that  would  yield  vertical  profiles  of  drop  size 
distributions  that  can  be  used  to  obtain  acceptable  estimates  of  vertical  profiles  of  ex¬ 
tinction,  backscatter,  and  absorption  in  very  low  stratus  clouds  and  their  associated  sub¬ 
cloud  regions.  (Very  low  clouds  are  considered  to  be  clouds  having  bases  less  than 
300  m  above  ground  level.  The  subcloud  region  is  the  region  between  ground  level  and 
the  base  of  the  cloud.)  In  addition  to  simulating  vertical  profiles  of  drop  size  distributions, 
the  model  also  produces  profiles  of  temperature,  pressure,  relative  humidity,  specific 
humidity,  mixing  ratio,  supersaturation,  liquid  water  content  (LWC),  potential  temperature, 
and  density.  There  will  be  a  user's  guide  (Kilmer,  1992)  for  this  model. 

The  report  is  arranged  as  follows.  First,  we  provide  a  brief  background  history  and  moti¬ 
vation  for  the  model.  This  is  followed  by  a  description  of  the  model  concept,  a  listing  of 
assumptions,  model  equations,  model  results  including  comparisons  with  field  data,  and 
a  summary  and  conclusions.  A  list  of  symbols  appears  at  the  beginning  of  the  text. 

2.  BACKGROUND 

In  the  late  1970's  the  Army  realized  that  the  vertical  characteristics  of  fog  and  very  low 
stratiform  clouds  were  not  sufficiently  well  characterized  for  EO  studies.  As  such,  field 
data  were  collected  during  periods  of  fog  and  low  stratus  clouds  in  Germany  and  England 
in  the  late  70's  and  early  80  s.  Each  set  of  data  consisted  of  vertical  profiles  of  drop  size 
and  drop  number  values  obtained  with  Knollenberg  FSSP-100  sensors  attached  to  bal¬ 
loons  that  ascended  and  descended  through  the  fogs  and  clouds.  Typical  LWC  profiles 
derived  from  these  data  are  shown  in  figure  1;  profiles  of  extinction  are  shown  in  figure 
2.  Note  that,  in  the  portion  of  these  figures  below  the  maximum  mass  loading  and  ex¬ 
tinction,  the  lower  parts  of  the  curves  are  concave  downward;  whereas  the  upper  parts 
are  convex.  Some  researchers  define  the  base  of  the  cloud  as  that  height  where  the 
inflection  of  log  (LWC)  versus  height  curves  occurs  as  shown  by  the  dashed  lines  in  figure 
1.  We  have  followed  this  definition  in  the  examp  es  modeled  for  this  report. 

The  first  set  of  data  collected  for  research  purposes  consisted  of  seven  profiles  of  drop 
size  and  drop  number  in  the  vicinity  of  Grafenwbhr,  Germany  (1976),  and  seven  profiles 
from  Meppen,  Germany  (1978).  These  daia  were  reported  by  Duncan  et  al.  (1980).  The 
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range  of  drop  radii  recorded  was  0.25  to  23.5  ^m.  A  second  set  (89  profiles  of  the  same 
size  range)  was  measured  at  Meppen  in  1980  as  reported  by  Lindberg  (1982).  A  third  set 
(5  profiles  of  drop  radii  0.2  to  50  ^m)  was  taken  at  Cardington,  England,  in  1983  (Lindberg 
et  al.,  1984).  The  data  from  Germany  were  generally  characteristic  of  very  low  stratus 
clouds;  whereas  those  from  England  were  obtained  from  higher  clouds. 


The  first  set  of  data  from  Germany  was  used  (Duncan  et  al.,  1980)  to  generate  the  basic 
empirical  model  as  follows.  The  logarithm  of  LWC  or  extinction  (oe)  for  each  of  three 
wavelengths  (0.55,  4.0,  and  10.6  /im)  at  each  height  z  +  20  m  was  plotted  against  the 
logarithm  of  LWC  or  oe  at  height  z,  producing  linear  relationships,  examples  of  which  are 
shown  in  figures  3  and  4.  A  mathematical  representation  of  figures  3  and  4  is 

x(z  +  20  m)  =  a'x(z)  +  b'  (1) 

where 

x(z  +  20  m)  =  In  LWC  or  ae  at  height  z  +  20  m 
x(z)  =  In  LWC  or  oe  at  z 

a'  and  b'  are  empirical  coefficients.  This  model  is  very  appealing  to  the  user  since  it  is 
simple,  very  easy  to  execute,  and  requires  only  a  single  value  of  LWC  or  oe  near  the 
surface  (2  m  above  ground)  as  input. 


Heaps  (1982)  and  Rachele  and  Kilmer  (1991),  in  addition,  showed  that  equation  (1)  can  be 
expressed  in  double  exponential  form,  that  is, 

LWC  or  =  A  exp  {  B  exp  Cz  }  (2a) 

where  z  is  expressed  in  meters.  A,  B,  and  C  are  related  to  a'  and  b'  from  equation  (1) 
according  to 


A 

B 

C 


(2b) 

(2c) 

(2d) 


and  yr  is  the  value  of  LWC  or  oe  at  a  reference  height  z  =  zr,  that  is, 

yr  ■  LWC  (zr)  or  oe  (zr)  (2e) 

Heaps’  also  found  that  the  Meppen  80  and  Cardington  cloud  data  yielded  profiles  that  are 
generally  consistent  with  Duncan's  earlier  results,  except  that  the  Cardington  data 
seemed  to  show  a  discontinuity  at  the  base  of  the  cloud.  Lindberg  and  Duncan  (1985) 
provided  adjustments  to  earlier  models  to  account  for  the  discontinuity  between  cloud 
and  subcloud  regions  for  the  higher  cloud  cases. 


The  LWC  and  extinction  models  of  Duncan  and  Heaps  are  basically  empirical  models  re¬ 
presenting  the  “average”  profiles  for  the  locales  where  the  data  were  collected.  In  con- 


1  verbal  conversation 
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trast,  we  are  hopeful  that  the  physical  model  developed  in  this  report  will  be  general 
enough  to  provide  “reasonable"  profiles  for  worldwide  applications,  while  at  the  same 
time  showing  why  one  could  expect  considerable  variation  in  drop  size,  extinction,  and 
LWC  profiles  on  a  worldwide  basis. 

3.  MODEL  CONCEPT  AND  ASSUMPTIONS 

Our  current  concept  of  this  model  envisions  an  ambient  atmosphere  that  is  characterized 
by  the  dynamics  of  its  mean  motion  as  modified  by  its  turbulent  elements,  but  has 
reached  a  steady  state  and  is  horizontally  homogeneous.  Imbedded  in  this  ambient  at¬ 
mosphere  is  a  layer  of  very  low  stratus  clouds.  We  assume  that  the  cloud  layer  is  capped 
by  an  inversion  layer.  In  addition,  we  also  envision  an  element  of  air  near  the  ground 
surface  having  a  density  sufficiently  less  than  the  ambient  air  to  rise  through  the  ambient 
air  up  to  the  top  of  the  cloud.  We  refer  to  this  element  of  air  as  a  representative  element. 
Other  elements  also  exist,  but  they  may  either  overshoot  or  undershoot  the  cloud  top 
height.  Of  all  the  elements  present,  we  focus  only  on  the  representative  element  in  this 
model  in  an  attempt  to  quantify  the  physical  and  thermodynamic  changes  that  take  place 
in  the  element  as  it  rises.  Imbedded  in  this  representative  element  of  air  is  a  population 
of  droplets.  We  refer  to  this  combination  of  air,  moisture,  and  droplets  as  a  droplet  clus¬ 
ter,  or  cluster  of  drops.  We  define  and  use  the  model  equations  in  a  way  that  character¬ 
izes  the  ambient  atmosphere  by  imposing  constraints  that  we  feel  are  reasonable.  The 
most  rigid  constraint  requires  that  the  model  equations  can  be  specified  by  using  a  small 
set  of  conventional  input  data  obtained  at  a  ground  level  reference  height  (say  2  m)-no 
profile  data  are  required. 

The  initial  drops  of  the  cluster  are  characterized  by  the  reference  level  drop  size  distrib¬ 
ution;  but  as  the  cluster  rises  the  drops  change  in  size,  and  therefore  the  distribution 
changes.  The  physical  and  thermodynamical  changes  of  the  drops  and  their  environment 
are  computed  as  though  a  hypothetical  completely  permeable  sphere  envelops  the  drops 
and  their  environment.  Hereafter,  the  terms  cluster  and  sphere  are  used  to  include  the 
drops  and  their  environment.  The  environment  outside  the  sphere  is  the  ambient  envi¬ 
ronment. 

The  model  allows  ambient  air  to  enter  the  sphere  as  it  ascends.  This  process  is  called 
entrainment  (Rogers,  1979;  Pruppacher  and  Klett,  1978,  1980).  Also  included  in  the  cur¬ 
rent  model  is  the  loss  of  water  (vapor  only)  from  the  sphere  by  a  process  we  call  shedding 
in  the  present  work.  This  is  consistent  with  our  notion  that  as  ambient  air  is  mixed  into 
the  sphere,  we  also  have  air  in  the  sphere  diffused  into  the  ambient  environment.  That 
is,  part  of  the  moisture  component  of  the  drop  environment  is  shed  into  the  ambient  air 
as  the  cluster  rises  through  the  ambient  environment.  Hence,  not  only  is  ambient  air 
being  entrained  into  the  cluster,  but  cluster  air  is  also  being  transferred  (shed)  simul¬ 
taneously  out  of  the  sphere  into  the  ambient  air. 

The  layer  from  the  base  of  the  cloud  to  the  user  input  value  for  the  top  of  the  cloud  is 
defined  as  the  cloud  layer.  The  base  of  the  cloud  is  defined  as  the  lowest  point  at  which 
the  cluster  environment  is  saturated  or  the  point  of  inflection  on  the  log  (LWC)  versus 
height  curve.  In  cases  where  there  are  such  a  large  number  of  such  hygroscopic  nuclei 
that  saturation  is  never  reached  or  is  delayed  considerably  beyond  what  would  reason¬ 
ably  be  considered  the  base  of  the  cloud,  the  latter  definition  is  useful. 
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The  microphysics  of  the  cluster  environment  is  computed  for  that  stage  when  the  cluster 
environment  is  still  distinguishable  from  the  ambient  environment.  After  an  initial  stage 
of  cloud  development,  the  cluster  and  ambient  environments  become  indistinguishable 
from  a  practical  measurement  standpoint.  At  this  stage  we  say  that  the  cloud  is  in  a  ma¬ 
ture  state,  and  the  cloud  and  subcloud  regions  are  in  steady  state. 

The  model  is  purposely  designed  to  require  a  minimum  number  of  reference  level  inputs 
(standard  observations)  and  a  knowledge  of  the  drop  size  distribution  at  the  reference 
height,  and  the  chemistry  of  the  dry  particles  upon  which  moisture  is  condensed.  We 
found,  however,  that  an  estimate  of  the  cloud  top  height  was  also  required. 

The  droplets  comprising  the  cluster  change  in  size  with  height  in  accordance  with  dif¬ 
fusion  controlled  differential  equations — the  smallest  remaining  in  approximate  equilib¬ 
rium  with  their  environment  and  the  largest  in  nonequilibrium.  Diffusion  of  water  vapor 
through  air  always  occurs  from  regions  of  higher  relative  humidity  to  regions  of  lower 
relative  humidity,  regardless  of  the  relationship  between  a  drop's  radius  and  its  critical 
supersaturation  value. 

The  ascent  rate  of  the  cluster  is  affected  by  buoyancy,  mass  of  liquid  water,  entrainment, 
and/or  forced  vertical  lifting.  The  cluster  stops  rising  in  the  vicinity  of  the  top  of  the  cloud. 

The  ambient  and  cluster  environments  have  the  same  relative  humidity.  This  assumption 
is  used  to  relate  the  water  vapor  mixing  ratio  in  the  ambient  environment  to  that  in  the 
cluster  sphere. 

The  temperature  difference  between  the  cluster  and  ambient  environments  decreases 
linearly  with  height  and  is  zero  at  the  top  of  the  cloud. 

The  ambient  environment  is  damp  and  nonadiabatic. 

The  drop  size  distribution  at  the  reference  height  is  defined  by  Shettle  and  Fenn  (1979) 
and  is  a  function  of  relative  humidity,  visibility,  and  geography  (urban,  rural,  maritime). 
However,  droplets  having  dry  radii  less  than  a  "lower  cutoff”  value  set  by  the  modeler  are 
omitted. 

The  chemistry  of  the  dry  particles  (condensation  nuclei)  is  adequately  characterized  by 
Hanel  and  Lehmann  (1981). 

The  equilibrium  size  of  a  drop  is  properly  defined  by  Hanel  and  Lehmann  (1981). 

The  cloud  and  subcloud  regions  are  connected,  and  hence  can  be  unified,  for  clouds 
whose  bases  are  less  than  300  m  above  the  ground. 

All  droplets  are  in  equilibrium  with  their  environment  inside  the  sphere  at  the  reference 
height. 

The  volume  of  the  sphere  that  envelops  the  droplets  is  equal  to  the  total  mass  of  dry  air 
inside  the  sphere  divided  by  the  density  of  this  dry  air.  The  total  masses  of  dry  air  and 
water  vapor  inside  the  sphere  are  affected  by  entrainment  and  shedding.  The  density  of 
dry  air  depends  on  its  partial  pressure  and  the  temperature  according  to  the  ideal  gas 
law. 
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The  total  pressure  of  damp  air  at  a  given  height  is  the  same  inside  and  outside  the 
sphere.  The  vertical  gradient  of  this  total  pressure  obeys  the  hydrostatic  approximation 
applied  to  damp  air  (dry  air  plus  water  vapor)  in  the  environment.  Since  the  water  vapor 
mixing  ratio  inside  the  sphere  may  be  different  from  the  ambient  environmental  value,  the 
partial  pressures  of  dry  air  and  water  vapor  inside  the  sphere  may  differ  from  the  corre¬ 
sponding  values  outside  the  sphere. 

The  concentration  of  drops  in  a  given  size  class  is  inversely  proportional  to  the  volume 
of  the  sphere. 

4.  FORMULATIONS 


4.1  Equilibrium  Drop  Size  Model  (Hanel) 

In  this  section  we  present  a  brief  summary  of  the  basic  formulations  that  relate  the  size 
of  a  droplet  to  the  relative  humidity  and  temperature  of  its  immediate  environment,  and 
the  chemical  composition  of  the  drop,  assuming  that  the  drop  is  in  thermodynamic  equi¬ 
librium  with  its  environment.  For  more  detail  one  should  read  H2nel  (1976),  who  initiates 
his  development  with  the  statement  that  his  model  is  based  on  an  equation  that  describes 
the  influence  of  curvature  and  dissolved  material  on  the  equilibrium  partial  vapor  pres¬ 
sure  Pr  over  an  uncharged  particle.  He  also  states  that  according  to  Harrison  (1965)  it  is 
necessary  to  relate  Pr  to  the  equilibrium  partial  vapor  pressure  Pw  over  a  plane  surface 
of  >  ;ier  being  saturated  with  dry  air.  Pw  is  calculated  for  the  total  pressure  P  of  the  moist 
air  surrounding  the  particle  and  for  the  particle's  absolute  temperature  T. 


After  some  discussion  Hdnel  concludes  that  the  ratio  of  Pr  to  Pw  can  be  written  as 


=  Aw  exp 


(3) 


where  Aw  is  the  water  activity  of  the  particle's  liquid  cover  with  a  hypothetical  plane  sur¬ 
face  in  absence  of  dry  air  in  the  gaseous  as  well  as  the  liquid  state,  and  can  be  measured 
within  a  closed  laboratory  system  at  well-defined  equilibrium  vapor  pressure  of  water  and 
temperature.  Also,  in  equation  (3),  o  is  the  surface  tension  on  the  particle's  surface,  vw 
the  specific  volume  of  pure  water  at  the  standard  pressure  of  1  bar,  r  the  radius  of  cur¬ 
vature  of  the  particle  surface,  and  Rv  is  the  specific  gas  constant  of  pure  water. 


The  exponent  part  of  equation  (3)  represents  the  curvature  correction  part  of  the  pressure 
ratio  when  enough  water  is  present  to  form  a  complete  shell  such  that  the  drop  can  be 
regarded  as  a  sphere.  This  curvature  correction  is  written  in  abbreviated  form  as 


Following  HSnel  and  Lehmann  (1981),  two  other  cases  also  are  considered  for  curvature 
correction:  (1)  a  particle  that  does  not  have  a  complete  liquid  cover  and  may  be  non- 
spherical  (expected  at  the  smallest  relative  humidities),  and  (2)  a  particle  in  a  transition 
region  between  the  two  extreme  cases.  The  three  cases  are  combined  in  the  following 
formulation: 
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(5) 


Assuming  that  the  particle  is  in  thermodynamic  equilibrium  with  the  surrounding  moist 
air,  Hanel  concludes  that  the  ratio  Pr/P w  of  equation  (3)  is  equivalent  to  the  relative  hu¬ 
midity  f  of  the  moist  air;  that  is, 

f  =  Pr/Pw  (6) 

and,  hence,  f  relates  to  Aw  and  S  as 

f  =  Aw<5  (7) 

which,  in  the  case  of  a  plane  surface,  reduces  to 

f  -  Aw  (8) 


At  this  point  it  is  appropriate  to  discuss  the  form  and  characteristics  of  Aw.  H3nel  states 
that  for  the  ideal  case  of  a  pure  solute,  that  is,  the  particle  is  an  aqueous  solution  of  a  pure 
solute,  Aw  can  be  written  in  four  different  equation  forms.  For  our  work  we  selected  the 
exponential  form 

Aw  =  exp  (  — (9) 


where  <j)  is  the  practical  osmotic  coefficient;  v  the  number  of  moles  of  ions  formed  from 
one  mole  of  solute,  being  equal  to  1  if  there  is  no  dissociation  of  the  solute  molecules; 
ns  the  number  of  moles  of  the  dry  solute;  and  nw  is  the  mole  number  of  water. 


Since  nw  =  mw/Mw  and  ns  =  ms/Ms,  where  mw  is  the  mass  of  water,  Mw  its  molar  mass, 
ms  the  mass  of  the  solute,  and  Ms  its  molar  mass,  equation  (9)  becomes 


Aw  = 


exp 


( 


—  (f)V 


Mw 

Ms 


(10) 


Defining  </>v 


Mw 

Ms 


Aw  =  exp 


as  r?s>  which  can  be  described  empirically,  equation  (10)  is  written  as 


(11) 


rf s  is  defined  as  the  exponential  mass  increase  coefficient. 


For  aerosol  particles  under  real  conditions  of  water  uptake,  HSnel  writes  equation  (11)  as 


Aw  — 


(12) 
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where 


mo  =  mu  +  ms  =  total  mass  of  a  dry  particle 
mu  =  mass  of  insoluble  material 
ms  =  mass  of  the  dry  mixed  solute 


The  exponential  mass  increase  coefficient  t)  is  evaluated  by  using  an  expression  that  is 
a  modification  by  Kilmer  suggested  by  Hanel  as  a  result  of  studies  by  Fitzgerald  et  al. 
(1982);  that  is. 


(13) 


where  the  polynomial  coefficients  aj  are  determined  by  curve  fitting  combined  data  from 
Hanel  and  Lehmann  (1981).  The  summation  in  braces  thus  represents  a  second  order 
effect  for  the  collection  of  all  measurements. 


The  combined  data  used  in  curve  fitting  include  two  categories  of  data  for  each  of  the  10 
samples  reported  by  Hanel  and  Lehmann  (1981):  (1)  discrete  data  for  Aw  <  0.965  and  (2) 
values  calculated  for  0.965  $AW^  1  from  Hanel  and  Lehmann's  approximate  extrapo¬ 
lation  formula;  that  is, 


(14) 


The  constants  ri°.  A,  B,  and  C  are  provided  by  Hanel  and  Lehmann  (1981)  for  each^of  their 
10  sets  of  measurements.  Since  t)°  —  the  limit  of  ij  as  Aw  1,  the  A,  B,  and  C  terms 
represent  the  difference  from  the  limiting  value. 

To  evaluate  the  curvature  part  of  equation  (3),  that  is,  equation  (4),  we  must  evaluate  the 
surface  tension  a,  T,  and  r.  The  expression  proposed  by  Hanel  for  evaluating  a  is 


a 


where 


(5) 

Two  a(T  To)  +  b  /  mw  \ 

) 

(Two  =  Uw(T0)  =  nw(T0  =  273.16  K)  =  75.6  dynes/cm 
a  =  0.153  dynes  cm"1  K"1,  -10°  C  <  T  <  30’  C 
b  =  equals  a  constant  (different  for  each  solute) 


(15) 


Hanel  (verbal  communication)  suggests  that  for  our  application  the  term  involving  b  is 
negligible. 


To  solve  for  r  we  assume  that  the  volume  of  a  droplet  consists  of  the  volume  of  the  dry 
particulate  plus  the  layer  of  water  condensed  onto  the  particle;  that  is, 


Taking  the  ratio  of  equations  (16a)  and  (16b)  gives 


(16a) 

(16b) 
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(17) 


or 


vw 

('3  -  ’?) 

vd 

rd3 

since  vw 
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mw 

Pw  (r3  -  '<?) 

md 

3 

Pdrd 

md 

Pdrd3 

mw 

Pw(r3  -  rd3) 

md 

“pT 


equation  (17)  can  be  expressed  as 


(18) 


(19) 


The  procedure  for  obtaining  T  will  be  discussed  in  later  sections  of  this  report. 

4.2  The  Vertical  Gradient  Expression  of  the  Equilibrium  Droplet  Relative  Humidity 

Of  necessity  for  later  developments  in  this  report  we  subscript  f  with  an  r  in  the  equations 
of  section  4.1.  In  effect,  if  a  droplet  is  in  equilibrium  with  its  immediate  environment  we 
designate  the  relative  humidity  as  fr.  Since  fr  was  shown  to  be  a  function  of  T  and  r,  for 
given  rd,  we  can  write  the  differential  of  fr  as 


d\{ 

dfr  -  dr  + 


JOl 

3T 


dT 


(20) 


If  we  then  assume  that  fr,  r  and  T  vary  with  height  z,  both  sides  of  equation  (20)  can  be 

multiplied  by  —  and  divided  by  dz  giving 
•r 


1  dfr 


d(  In  fr) 


fr  dz  dr 

For  convenience,  we  write 
d  In  fr 


dr 

dz 


+ 


3(lnfr)  qj 


dT 


dz 
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dr 


—  d  In  fr 

_  Q  and  - — —  =  H. 

o  T 


4.3  Nonequilibrium  Model  (Rachele  and  Kilmer) 

We  consider  two  general  properties  affecting  all  droplets:  the  change  in  mixing  ratio  and 
thermodynamics  of  the  cluster  environment  as  it  ascends.  The  following  symbols  are 
used  in  this  section: 

ote  =  entrainment  parameter 

R  =  radius  of  sphere 

Wv  =  grams  of  water  vapor  per  gram  of  dry  air  in  sphere 
W'v  =  grams  of  water  vapor  per  gram  of  dry  air  in  ambient  air 
z  =  height  above  ground  level 

b  =  an  adjustable  parameter 

c  =  an  adjustable  parameter 

e  =  ratio  of  molecular  weight  (MW)  of  water  to  mean  MW  of  dry  air  =  0.621972 


22 


V  =  temperature  of  the  ambient  air  in  kelvins 
T  =  sphere  temperature  in  kelvins 

Wl  =  mass  loading  (liquid  water  plus  dry  material)  mixing  ratio  inside  the  sphere 
Zt  =  height  above  ground  at  top  of  cloud 

Le  —  specific  enthalpy  of  vaporization  of  water  at  temperature  of  sphere 

a  *=  temperature  radiation  constant 

Cp  =  specific  heat  of  dry  air  at  constant  pressure 

CpG  =  specific  heat  of  water  vapor  at  constant  pressure 

S[_  =  specific  heat  of  liquid  water 

Ssj  =  specific  heat  of  solid  or  dissolved  substance  i 

Wsi  =  grams  of  solid  or  dissolved  substance  i  per  gram  of  dry  air 


The  first  of  two  basic  equations  used  in  this  development  (discussed  in  Pruppacher  and 
Klett  (1978)  for  saturated  conditions)  describes,  in  differential  form,  the  thermodynamics 
of  an  ascending  “parcel.”  This  equation  has  been  modified  in  the  present  work  and  ex¬ 
tended  to  our  concept  of  a  rising  cluster  (see  section  4.4)  to  become 

Le  dWv  dT  _  gT(1  +  Wv)  f  1  +  WJc  } 

C*  dz  dz  CpT'  (  1  +  Wv/c  j  (22a) 

+  S  +  E-|-Et  +  Ac-|-  (Turb) 


where  E  is  the  "horizontal"  entrainment, 


E  = 


Le(Wv-W') 

c; 


T') 


r  3«eT  ] 

4RT' 


~  1  -I-  W  Jc  ~ 
_  1  +  W  Jr.  _ 


(22b) 


The  entrainment  parameter  ae  is  discussed  in  section  4.4. 


S  is  shedding:  This  effect  may  be  considered  analogous  to  “horizontal"  entrainment  re¬ 
presented  by  equation  (22b)  (but  with  mass  flow  in  the  opposite  direction).  However,  the 
Wv  and  T  of  the  material  being  shed  equal  the  average  values  in  the  sphere.  Hence,  since 
this  would  be  analogous  to  having  Wv  =  Wv'  and  T  =  T'  in  equation  (22b),  S  =  0. 


Ej  is  the  "vertical"  entrainment  from  inversion  layer  above  the  cloud  top, 

-C(ZT  -7)* 


Et  - 


be 


(22c) 


Cp 

Ac  is  the  net  radiation  effect  in  the  cloud  and  is  an  unknown  function  of  z, 

Ac  -  F(z)  (22d) 

(Turb)  is  the  turbulence  effect,  after  Priestley  (1953), 

(Turb)„  (22e, 

E  and  (Turb)  are  not  used  simultaneously  in  equation  (22a)  since  they  represent  the  same 
effect  but  are  viewed  from  different  perspectives. 
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Cp  is  the  constant  pressure  heat  capacity  of  a  mass  that  is  in  the  cluster  sphere  and  that 
includes  exactly  1  g  of  dry  air.  Cp  is  the  analogous  quantity  for  a  mass  that  is  in  the  am- 


bient  environment. 

c;  -  Cp  +  Wv  4-  SLWL  +  ISsiWsl 

(22f) 

Cp'  =  Cp  +  CpoW^ 

(22g) 

In  simulations  using  the  above  equations,  "horizontal”  entrainment  or  the  Priestley  (1953) 
turbulence  effect  tends  to  increase  the  magnitude  of  dT/dz  when  dT/dz  is  negative,  but 
“vertical”  entrainment  from  a  warmer  layer  above  the  cloud  decreases  the  magnitude  of 
dT/dz  while  this  gradient  is  negative  and  eventually  forces  dT/dz  to  become  positive. 
Hence,  "horizontal”  entrainment  and  turbulence  tend  to  increase  LWC,  unless  the  amount 
of  liquid  water  that  must  be  evaporated  to  bring  the  mixing  ratio  of  drier  entrained  air  up 
to  that  of  damp  air  in  the  sphere  exceeds  the  amount  of  liquid  water  gained  by  the  effect 
of  entrainment  on  dT/dz;  however,  “vertical”  entrainment  from  an  inversion  layer  de¬ 
creases  LWC. 


Different  amounts  of  entrainment  can  be  modeled  by  adjusting  ae  in  equation  (22b).  We 
are  not  aware  of  a  theoretical  value  or  a  theoretical  expression  that  would  generate  a 
value  for  ae-  For  the  very  low  stratus  clouds  that  we  have  studied  to  date,  ae  =  0.03 
appears  to  be  a  reasonable  value.  Our  current  equation  for  approximating  Et  is  given  as 
equation  (22c).  Equation  (22e)  is  linear  in  the  difference  T  -  T',  which  in  turn  is  modeled 
as  being  linear  with  z  -  zr,  the  height  above  the  reference  level.  It  is  assumed  that 
equation  (22d)  can  be  represented  by  an  equation  of  the  form 

Ac  =  a(z  -  zr)  (22h) 


Much  of  the  justification  for  including  the  form  given  as  the  right-hand  side  of  equation 
(22h)  is  the  empirical  success  enjoyed  in  matching  observed  LWC  profiles  with  profiles 
simulated  when  this  term  is  included.  Recently  Klett*  found  this  term  to  be  a  reasonable 
approximation  of  gross  effects  of  radiation.  As  a  result  of  the  above,  equation  (22a),  in 
terms  of  entrainment,  becomes 


Le  dWv  dT 
r*  dz  dz 


+wj  r  l 

* 

'  T 

' P 


'4w,-w,)+£l(t-t') 


1  -I-  Wv/c 
1  +  W  Jr. 


-c(zT  -zy 


+  a(z-zr) 


(In  this  study,  the  LSs,WSi  term  is  set  equal  to  zero  because  the  contribution  from  this 
relatively  small  amdunt  of  matter  is  negligible.) 


The  second  basic  equation  describes  the  changes  in  the  mixing  ratio  of  an  ascending 
“parcel.” 


*  Klett,  James  D  .  unpublished  report 
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Pruppacher  and  Klett  (1978)  express  the  mixing  ratio  (Wv)  in  terms  of  supersaturation 
(Sv,w).  saturation  vapor  pressure  (es),  and  partiaf  pressure  of  dry  air  (P<j)  as 

W  v  =  (1+Sv,w)€-^_  (23) 

where 

Sv.w  -  f-1. 


After  equation  (23)  is  differentiated,  some  quantities  are  replaced  by  using  the  Clausius- 
Clapeyron  expression  (while  neglecting  the  volume  of  liquid  water  produced  by 
condensation,  as  neglecting  this  volume  causes  only  about  5  parts  per  million  error  at  a 
temperature  of  273.16  K) 


des  =  lees  dT 
dz  RvT2  dz 


(24) 


and  the  hydrostatic  approximation  applied  to  the  total  pressure  in  the  environment  (see 
derivation  in  the  appendix) 


dPd  _  ~0pd 
dz  "  RdT' 


dWv 

dz 


(25) 


Equation  (23)  in  differential  form,  in  terms  of  dSv,w/dz,  becomes 

= r,  +  s  iT-! _ +  s  iT  l’dT  i  °  (  1+w»  V 

dz  [  »-*]  L  wv  (.  +  W„)  J  dz  [  +  R>Tzd2  +  RdT'  (  ,  +  w;/t  J 

Equation  (26)  represents  a  significant  part  of  droplet  growth  as  described  below. 


(26) 


The  particles  in  nonequilibrium  grow  in  accordance  with  Maxwell's  diffusion  equation, 
which  may  be  written  as 


—  =  —  (f  _  f ) 

dt  r  ^  'rf 

(27a) 

or,  as  we  prefer 

=  (f  _  f ) 
dz  rw  "  'rf 

(27b) 

where  w  is  the  cluster  ascent  rate,  and 

f  Pw*-e  /  i-e  \") 

-  \  Dvd,  +  KT  (  R„T  ’j) 

(27c) 

To  evaluate  equation  (27b),  one  must  have  estimates  of  f  and  fr  as  functions  of  z.  We 
approximate  them  by  using 


25 


(28a) 


fr  =  fr1  +  -Jj“  AZ 

where  the  subscript  (1)  relates  to  a  height  immediately  below  z,  that  is,  z  -  Az. 
From  equation  (26),  we  obtain 

df  _f[~  1  dWv  l-e  dT  g  (  1  +  Wv  \ 

dz  Y  Wv(1  4-  Wv/c)  dz  RJ2  dz  RdT'  ^  1  +  w'v/£  J 


(28b) 


(29) 


and  from  equation  (21),  we  obtain 


[g-5T-h 


Substituting  equations  (28)  to  (30)  into  equation  (27b)  gives 

dr  DAz  I"  (fi  ~  fn)  f  dWv  (. .  ,  M  \  _dT 

dz  rw  Az  +Wv(1+Wv/r.)  dz  +  dz 


Solving  equation  (31)  for  dr/dz  gives 

_dL  =  Y  f  (f1~fn)  ,  f  dWv  ,  /,,,  LJ  \  dT  __gf_  /  1+w'v  \ 

dz  Az  Wv(1+wv/r.)  dz  RyT2  ;  dz  RdT'\^1  +  w J 

where 


(30) 


(31) 


(32) 


(33) 


Examination  of  equation  (32)  shows  that  dr/dz  is  governed  by  equilibrium  factors  for  small 
drops  since  G  is  approximately  3000  for  r  =  3  x  10  "5  cm;  hence,  DAzfrG/(rw)  is  approxi¬ 
mately  200  for  Az  =  10  cm,  D  =  6.16  x  10“7,  and  w  =  3  cm/s.  In  contrast,  when  r  is  1  x 
10~3  cm,  G  is  approximately  75;  hence,  DAzfrG/(rw)  is  approximately  0.15.  Comparison 
of  the  latter  value  with  1  shows  that  growth  of  the  larger  drops  is  controlled  largely  by  the 
rate  of  diffusion  of  water  vapor.  Therefore,  the  larger  drops  are  not  expected  to  be  at  their 
equilibrium  sizes. 


We  initiate  the  solution  of  equation  (32)  by  assuming  that  all  drops  are  in  equilibrium  at 
z  =  2  m,  so  that  f  -  fr  =  0  for  all  drops.  Hence,  at  one  step  above  the  reference  height 
(2  m  +  Az),  dr/dz  is  evaluated  for  fi  -  fr1  =  0  for  all  drops.  Above  the  reference  height. 
Az  is  generally  set  equal  to  10  cm. 
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For  each  drop  size  class,  r(z  +  Az),  the  drop  radius  at  height  z  +  Az,  typically  is  calcu¬ 
lated  by 


r(z  +  Az)  =  r(z)  + 


Az  /  dr 
2  ^  dz 


A.  ) 

dZ  r«z+Az ) 


(34) 


during  the  final  predictor-corrector  iteration  for  a  given  height  step.  However,  during 
rapid  evaporation,  there  is  some  danger  that  the  last  term  on  the  right  side  would  have 
a  sufficiently  large  negative  value  to  compromise  the  accuracy  of  the  calculation  and  even 
lead  to  numerical  instability.  Therefore,  no  drop  radius  is  ever  allowed  to  decrease  by 
more  than  half  of  the  difference  between  that  drop  radius  (at  height  z)  and  the  dry  particle 
radius  for  the  same  size  class  during  a  single  height  step  (Az). 


4.4  Cluster  Ascent  Rate 

The  cluster  ascent  rate  w  is  approximated  using  formulations  that  consider  the  effects  of 
buoyancy,  mass  loading,  entrainment,  and  turbulence.  The  Coriolis  force  is  not  included. 
The  general  form  of  the  cluster  ascent  equation  is 

w  -  g  ( -  l)  -  E  -  k,w  (35) 

where  the  first  term  on  the  right  includes  the  effects  of  buoyancy  and  mass  loading,  the 
second  is  due  to  entrainment,  and  the  third  term  is  due  to  turbulence,  in  line  with 
Priestley.  As  was  done  in  section  4.3,  we  allow  either  entrainment  or  Priestley  turbulence 
to  be  activated  in  equation  (35),  but  not  both.  If  entrainment  is  activated  in  equation  (35), 

-kiw  would  be  omitted,  and  E  would  be  calculated  as  ■—  w2.  R  is  the  radius  of 

the  sphere,  the  initial  value  of  which  is  fixed  at  2  m  since  we  generally  assume  the  ref¬ 
erence  height  to  be  2  m  above  ground.  The  radius  of  the  sphere  is  affected  only  by  the 
change  of  density  of  dry  air  in  the  sphere. 

In  obtaining  an  expression  for  the  density  ratio  p'/p  in  equation  (35),  the  following  items 
are  used  from  assumptions  8  and  17  in  section  3:  It  is  assumed  that  the  total  air  pressure 
of  the  cluster  environment  is  equal  to  the  total  pressure  of  the  ambient  air.  It  also  is  as¬ 
sumed  that  the  ambient  air  contains  no  liquid  water  or  solid  matter. 


Since  hydrostatic  equilibrium  is  assumed  to  hold  for  the  total  pressure  and  density  of 
ambient  air,  the  ratio  of  the  density  p'  of  ambient  air  to  the  mean  density  p  of  the  cluster 
environment  is  given  by 

P'  T(H-W</t)(l  +  W») 

p  T'(1  +  W,  +  We)  (1  +  W,/r) 

Our  formulation  of  entrainment  (which  is  represented  by  the  E  term  in  equation  (35)) 
considers  the  entrainment  of  ambient  air  into  the  rising  cluster  only;  it  does  not  account 
for  entrainment  of  dry  stable  air  above  the  cloud  into  the  cloud.  (However,  the  latter  effect 
is  represented  by  the  last  term  in  the  cluster  lapse  rate  equation.)  The  formulation 
evolves  as  follows:  A  rising  cluster  (including  intimately  associated  damp  air)  displaces 
ambient  air.  The  displaced  air  from  above  the  cluster  fills  the  void  below,  except  that  a 
small  fraction  ae  of  the  displaced  mass  is  “entrained"  into  the  cluster.  The  volume  of 
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displaced  air  is  equal  to  rcR*dz,  where  R  =  radius  of  the  hypothetical  sphere  that  defines 
the  cluster  environment  and  dz  is  an  infinitesimal  increment  of  height.  Therefore,  the 
volume  of  air  entrained  while  the  cluster  ascends  by  this  height  increment  dz  is 
aertR*dz.  Multiplying  this  volume  by  p\  the  density  of  damp  air  in  the  environment,  and 
dividing  by  the  mass  of  the  cluster  (including  intimately  associated  damp  air)  47tR3p/3 
shows  that 


1  dm  3  «e  P' 
m  dz  ARP 


(37) 


The  work  required  to  accelerate  a  given  amount  of  mass  Mi  from  a  velocity  of  (u,v,0)  to 
(u,v,w)  is  taken  to  be  equal  to  MiW*/2,  the  gain  in  kinetic  energy.  If  the  change  in  height 
required  for  this  increase  in  velocity  is  h,  and  the  force  required  for  this  acceleration  is 
assumed  to  be  constant,  then  the  force  is  MiW*/(2h).  However,  the  total  mass  entrained 
during  this  change  in  height  is  h(dm/dz).  Substituting  this  expression  for  Mi  shows  that 
the  total  force  is  (w*/2)dm/dz.  Dividing  by  the  total  mass  of  the  cluster  and  then  substi¬ 
tuting  for  (1/m)dm/dz  shows  that  the  force  necessary  to  accelerate  the  mass  being  en¬ 
trained  per  unit  of  total  mass  in  the  cluster  is 


c  1  dm  w?  3  ae  ft'  2 
m  dz  2“8RpW 


(38) 


This  is  included  as  part  of  the  E  term  in  equation  (35).  One  of  the  equations  used  in 
solving  for  vertical  gradients  of  temperature  and  water  vapor  mixing  ratio  in  the  cluster 
is  obtained  by  extending  a  first  law  of  thermodynamics  equation  from  Pruppacher  and 
Klett  (1980)  to  include  unsaturated  conditions  and  condensed  matter  in  the  cluster  and 
our  approximation  of  a  turbulence  effect.  When  the  entrainment  expression  from  the 
present  work  is  included,  the  resulting  equation  is 


<WV  dT 
r’  dz  dz 


gT(i  +  w„) 


C  T' 


L.(WV  -  W„)  Cp 

- —  +  ~  (T  -  T') 


(39) 


1  +  Wv/c 
1  +  W ' /t 


be 


-C(ZT  -Z>* 


+  a(z  -  zr) 


4.5  Temperature  and  Potential  Temperature 


The  temperature  of  the  cluster  environment  varies  in  accordance  with  the  differential  form 
given  by  equation  (39).  Note  that  it  is  also  a  function  of  the  mixing  ratio  of  the  cluster 
environment  and  the  entrainment  of  ambient  air.  Based  on  our  experience  and  many 
trials,  we  assume  that  the  ambient  temperature  for  a  mature  cloud  phase  relates  to  the 
cluster  temperature  as  follows 


This  equation  states  that  the  temperature  difference  between  the  cluster  and  ambient 
environments  varies  linearly  with  height,  and  equals  zero  at  the  top  of  the  cloud. 

The  potential  temperature  of  the  cluster  is  stated  in  the  traditional  form  as  presented  in 
R.  R.  Rogers  (1979,  1989);  that  is. 
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dfl 

dz 


Le  dWv 
o*  dz 


(Wv  -  W*w)  +  (T  -  T') 


1  dm 
m  dz 


(41) 


The  potential  temperature  of  the  ambient  environment  is  based  on  the  expression  of  the 
potential  temperature  of  damp  air;  that  is, 

(p  \ 

-jr)  (42) 

where  Rm  =  R<j(1  +  0.61q') 

Cpm  =  Cpd(1  +  0.84q') 

Pc  =  1000  millibars  or  1  x  106  dynes  cm-2 

In  differential  form  equation  (42)  may  be  approximated  by 


T  d0' 
0'  dz 


dr 

dz 


■'pm 


+  r  !«(£)(■ 


0.61Rd 


'pm 


0-84CptJRrn 


'pm 


dq' 

dz 


(43) 


4.6  Bulk  Phase  Change  Relationship 


For  inclusion  in  a  matrix  equation  that  is  appropriate  at  the  reference  height  only,  the 
phase  change  of  moisture  from  vapor  to  liquid  water  is  based  on  the  change  of  mixing 
ratio  of  the  cluster  air  assuming  that  total  moisture  is  conserved;  that  is,  dWL  =  -dWv. 
(This  simple  equality  is  not  appropriate  when  entrainment  and  shedding  are  included.) 
To  maintain  consistency  with  other  equations,  we  use  an  approximate  form  of  equation 
(26)  expressed  in  terms  of  f;  that  is. 


1  df  1  dWv 


Le  dT 


f  dz  W„  dz 


RVT 


2  dz 


RdT 


(44) 


Above  the  reference  height,  f  is  calculated  directly  as 
W„  P„ 


f  = 


r.  e. 


(45) 


where  es  is  calculated  as  a  function  of  temperature  according  to  Lowe  and  Ficke  (1974). 


4.7  Mass  (Moisture)  Conservation  (Kilmer  and  Rachele) 

The  following  assumptions  are  made: 

a.  The  volume  of  the  sphere  defined  by  the  cluster  equals  the  total  mass  in  the  cluster 
sphere  divided  by  the  overall  density  of  the  cluster  sphere. 

b.  The  number  of  particles  (droplets)  N,  per  cubic  centimeter  is  inversely  proportional 
to  the  volume  of  the  cluster  sphere. 

c.  This  mass  balance  equation  for  total  water  must  be  satisfied:  The  sum  of  (1)  the  mass 
of  water  vapor  in  the  sphere  at  height  z  plus  (2)  the  mass  of  liquid  water  in  the  sphere 
at  height  z  plus  (3)  the  mass  of  water  vapor  gained  by  entrainment  during  height  step 
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dz  must  equal  the  sum  of  (4)  the  mass  of  water  vapor  in  the  sphere  at  height  z  +  dz 
plus  (5)  the  mass  of  liquid  water  in  the  sphere  at  height  z  +  dz. 

d.  The  rate  of  change  with  height  of  total  water  mixing  ratio  in  the  sphere  is  equal  to  the 
rate  of  change  of  the  vapor  mixing  ratio  plus  rate  of  change  of  mixing  ratio  of  liquid 
water;  that  is, 

_dj_  =  dW^  dWL 

dz  dz  dz 

To  help  in  expressing  these  assumptions  in  the  form  of  an  equation,  a  variable  *  is  de¬ 
fined  as 


X 


(46) 


The  total  mass  of  dry  air  in  the  sphere  is 
A  3 

ma  Q  ^  ^  Pa 


(47) 


Solving  equation  (46)  for  R  and  substituting  into  equation  (47)  gives 

4  3 

ma  j  n  X  Par 


(48) 


Taking  the  natural  logarithm  of  both  sides  and  then  taking  derivatives  with  respect  to  z 
yields 


1  dma  _  3  dx 
ma  dz  X  dz 


f49) 


d* 

This  equation  shows  (by  multiplying  both  sides  by  x3)  that  3x2-^-  may  &e  replaced  with 
X 3  dma 
ma  dZ 

Using  the  definition  of  x.  assumption  3  is  represented  by  the  following  equation: 

N„ 


y  *  X3ParWv  +  (  y  *)  pw  R?  ^  Nir  (rf  -  4)  +  47t  X2ParWvd*  = 

2  Nte 

y^(x  +  dx)3Par(Wv  +  dWv)  +  (-*-*)  PwRfXNir[(ri  +  dr.)3-ry  -  T  *RVadl 

i=1 

Starting  with  this  equation  and  then 


(50) 


a.  dividing  through  by  y  n  dz,  discarding  the  second  and  third  powers  of  differentials, 
and  collecting  terms,  J 
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b.  replacing  3*2  with  ^7  ^ 


c.  replacing  x 3  with 


R3Pa 

Par 


,  ,  .  1  dma  .iL.  3aePa' 

d.  replacing  — - r—  with 


ma  dz 


4Rpa 


and 


e.  multiplying  through  by  —  n  and  rearranging  the  resulting  equation  (to  facilitate  identi- 

O 

tying  individual  terms  with  physical  quantities)  produces 

N„ 


4  d3 
-  y  *R  Pa 


dWv 

dz 


2  '  x  (4«) 

-  «R2«ePa  *(Ww  -  Ww)  +  R 


Iml 


(51) 


This  equation  may  be  interpreted  as:  The  mass  of  water  made  available  by  decrease  in 
water  vapor  mixing  ratio  with  height  in  the  sphere  equals  the  sum  of  masses  of  water 
used  to  increase  the  water  vapor  mixing  ratio  of  entrained  air  to  that  in  the  sphere  and  to 
increase  the  amount  of  liquid  water  in  the  sphere. 


Multiplying  both  sides  of  this  equation  by 


-3 

4rrrt3 


and  rearranging  again  produces 


/  Rr  \3  2  dr,  dWv  3a-  . 

«»(  x  dT  +  »•  ST  ~  "  p‘  (w«  “  w»' 

'  '  1=1 


(52) 


The  mass  balance  equation  in  this  form  is  represented  by  the  next-to-last  row  of  the  ma¬ 
trix  equation  that  is  solved  at  each  height  level. 


4.8  Drop  Size  Distribution 


We  assume  that  the  drop  size  distribution  at  the  reference  height  (2  m)  is  adequately  de¬ 
fined  by  the  Shettle-Fenn  (1979)  formulation  as  discussed  in  this  section. 


Shettle  and  Fenn  use  a  bimodal  lognormal  distribution 


n(r) 


dN(r) 

dr 


where 


I 

i=1 


N, 

r<r((27r)1/2  ln(10) 


( log  r  -  log  r,)2 


N(r)  =  the  cumulative  number  density  of  particles  of  radius  r 
oj  =  the  standard  deviation 

rj  =  the  mode  radius 

Nj  =  the  number  density  corresponding  with  rj 


(53) 


Four  different  aerosol  models  for  the  atmospheric  boundary  layer  near  the  earth's  surface 
were  developed  by  Shettle  and  Fenn.  The  models  differ  in  particle  size  distribution  and 
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particle  refractive  index.  Table  1  lists  the  parameters  defining  the  size  distributions  in 
accordance  with  equation  (53)  for  these  models. 

The  choices  of  N  in  table  1  were  normalized  by  Shettle  and  Fenn  to  correspond  to  1 
particle/cm3.  The  actual  size  distributions  can  be  renormalized  for  the  visibility  being 
used,  or  one  can  use  a  distribution  from  those  characterized  by  tables  3,  4,  and  5  for  rural, 
urban,  and  maritime  conditions,  respectively. 

TABLE  1.  CHARACTERISTICS  OF  THE  AEROSOL  MODELS  OF  THE 
LOWER  ATMOSPHERE  (SHETTLE  AND  FENN,  1979) 


Aerosol  Model 

Size  Distribution 

Type 

Nj 

n* 

°i 

RURAL 

0.999875 

003 

0.35 

Mixture  of  water-soluble 

0.000125 

0.5 

0.4 

and  dust-like  aerosols 

URBAN 

0.999875 

0.03 

0.35 

Rural  aerosol  mixture 

0.000125 

0.5 

0.4 

with  soot-like  aerosols 

MARITIME 

Continental 

Origin 

1. 

0.03 

0.35 

Rural  aerosol  mixture 

Oceanic 

1. 

0.3 

0.4 

Sea  salt  solution 

Origin 

in  water 

TROPOSPHERIC 

1. 

0.03 

0.35 

Rural  aerosol  mixture 

*  These  mode  radii  correspond  to  moderate  humidities  (70  to  80  percent);  values  of  rj  as 
function  of  humidity  are  given  in  table  2. 
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TABLE  2.  MODE  RADII  FOR  THE  AEROSOL  MODELS  AS  A  FUNCTION  OF 
RELATIVE  HUMIDITY  (SHETTLE  AND  FENN,  1979) 


Relative 

Humidity 

(%) 

Tropospheric 

ri 

Rural 

ri  r2 

Maritime 

ra 

Urban 

ri  r2 

0 

0.02700 

0.02700 

0.4300 

0.1600 

0.02500 

0.4000 

50 

0.02748 

0.02748 

0.4377 

0.1711 

0.02563 

0.4113 

70 

0.02846 

0.02846 

0.4571 

0.2041 

0.0291 1 

0.4777 

80 

0.03274 

0.03274 

0.5477 

0.3180 

0.03514 

0.5805 

90 

0.03884 

0.03884 

0.6462 

0.3803 

0.04187 

0.7061 

95 

0.04238 

0.04238 

0.7078 

0.4606 

0.04904 

0.8634 

98 

0.04751 

0.04751 

0.9728 

0.6024 

0.05996 

1.1691 

99 

0.05215 

0.05215 

1.1755 

0.7505 

0.06847 

1.4858 

The  radii  in  the  above  table  are  in  micrometer  {/xm)  units. 


TABLE  3.  TOTAL  NUMBER  DENSITY  FOR  THE  RURAL  AEROSOL  MODEL  AS  A  FUNCTION 

OF  RELATIVE  HUMIDITY  AND  METEOROLOGICAL  RANGE  (SHETTLE  AND  FENN,  1979) 


Relative 


Humidity 

~ 

0% 

50% 

70% 

80% 

90% 

95% 

98% 

99% 

Visibility 

N1 

= 

2  Km 
19953058 

192537.16 

178656  10 

129883  79 

85526.28 

68260.42 

4867499 

37244.91 

N2 

24.94 

24.07 

22  33 

16  24 

10.69 

8.53 

6.09 

4.66 

Visibility 

N1 

= 

5  km 
79076.33 

76304  75 

70803.53 

51474  48 

33895  07 

2705241 

19290.47 

14760  60 

N2 

= 

9.89 

9  54 

8.85 

6  44 

4.24 

3.38 

2.41 

1.85 

Visibility 

N1 

: 

10  km 
38924  91 

3756062 

34852  67 

25338  05 

16684.67 

13316.41 

9495  64 

7265.83 

N2 

= 

4.87 

4  70 

4.36 

3.17 

209 

1.66 

1.19 

0.91 

Visibility 

N1 

: 

23  km 
16230.63 

15661.76 

14532.62 

10565.23 

6957.06 

5552  58 

3959.42 

302965 

N2 

= 

2.03 

1.96 

1.82 

1.32 

0.87 

0.69 

049 

0.38 

Visibility 

N1 

= 

50  km 
6803  78 

6565.31 

6091.98 

4428  90 

2916.35 

232761 

1659.76 

127601 

N2 

= 

0.85 

0.82 

0.76 

0.55 

0.36 

0.29 

0.21 

0.16 

TABLE  4.  TOTAL  NUMBER  DENSITY  FOR  THE  URBAN  AEROSOL  MODEL  AS  A  FUNCTION 

OF  RELATIVE  HUMIDITY  AND  METEOROLOGICAL  RANGE  (SHETTLE  AND  FENN,  1979) 


Relative 

Humidity 

= 

0% 

50% 

70% 

80% 

90% 

95% 

98% 

99% 

Visibility 

N1 

= 

2  km 
220039.51 

210325.47 

163580.19 

107935.91 

69877.38 

45967.45 

26535.53 

18383.38 

N2 

— 

27.51 

26.29 

20.45 

13.49 

8.74 

5.75 

3.32 

2.30 

Visibility 

N1 

= 

5  km 
87204.26 

83354.47 

64828.76 

42776.28 

27693.23 

18217.44 

10516.34 

728535 

N2 

= 

10.90 

10.42 

8.10 

5.35 

3.46 

2.28 

1.31 

0.91 

Visibility 

N1 

= 

10  km 
42925.84 

41030.80 

31911.62 

2105640 

13631.85 

8967.44 

5176.62 

3586.27 

N2 

~ 

5.37 

5.13 

3.99 

2.63 

1.70 

1.12 

0.65 

0.45 

Visibility 

N1 

= 

23  km 
17898.91 

17108.73 

13306.28 

8779.95 

5684.11 

3739  18 

2158.51 

1495.38 

N2 

= 

2.24 

2.14 

1.66 

1.10 

0.71 

047 

0.27 

0.19 

Visibility 

N1 

; 

50  km 
7503.11 

7171.87 

5577  91 

3680.50 

2382.74 

1567.44 

904.83 

62685 

N2 

= 

0.94 

0.90 

0.70 

0.46 

0.30 

0.20 

0.11 

008 

TABLE  5.  TOTAL  NUMBER  DENSITY  FOR  THE  MARITIME  AEROSOL  MODEL  AS  A  FUNCTION 

OF  RELATIVE  HUMIDITY  AND  METEOROLOGICAL  RANGE  (SHETTLE  AND  FENN,  1979) 


Relative 

Humidity 

= 

0% 

50% 

Visibility 

N1 

_ 

2  km 
96516.79 

88640.83 

N2 

= 

974  92 

895.36 

Visibility 

N1 

; 

5  km 
3825075 

35129.41 

N2 

= 

386.37 

354.84 

Visibility 

N1 

: 

10  km 
18828.73 

17292.27 

N2 

= 

190.19 

174.67 

Visibility 

N1 

: 

23  km 

7851  07 

7210.41 

N2 

= 

79.30 

72.83 

Visibility 

N1 

= 

50  km 
3291.12 

3022  56 

N2 

= 

33.24 

3053 

70% 

80% 

90% 

70038.09 
707  46 

35078.01 

354.32 

24467.13 

247.14 

27756.92 

280.37 

13901  83 
14042 

9696.61 

97.95 

13663.20 

138.01 

6843  10 
69  12 

4773.11 

48.21 

5697.18 

57.55 

2853.39 

28.82 

1990.26 
20  10 

2388  22 
24.12 

1196.12 

1208 

834.30 

8.43 

95% 

98% 

99% 

17603.12 

177.81 

11000.60 

111.12 

7438.74 

75.14 

6976.32 

70.47 

4359.67 

44.04 

2946.06 

29.78 

3434.06 

34.69 

2146.02 

21.68 

1451.17 

1466 

1431.91 

14.46 

894.83 

9.04 

605.10 

6.11 

60025 

6.06 

375.11 

3.79 

253.65 

2.56 
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The  Shettle-Fenn  distribution  has  three  constraints;  that  is,  the  distribution  depends  on 
visibility,  relative  humidity,  and  aerosol  origin  (maritime,  urban,  or  rural). 

The  selected  distribution  is  used  to  determine  the  number  of  droplets  in  each  size  class 
as  defined  at  the  reference  level.  Lower  and  upper  bounds,  ra  and  rt>,  respectively,  are 
set  for  the  droplet  radius  in  each  size  class.  The  number  of  droplets  in  a  given  size  class 
is 


where  n(r)  is  defined  in  equation  (53)  at  the  reference  level.  Above  the  reference  level, 
the  number  of  droplets  per  cubic  centimeter  in  a  given  size  class  is  inversely  proportional 
to  the  volume  of  th?  sphere.  At  the  reference  level,  the  representative  radius  for  that  size 
class  is  calculated  as  the  mean  volume  radius  for  ra<,r  <,  rb;  that  is, 


n(r)  dr  / 


1/3 


This  form  was  chosen  because  of  the  importance  of  r3  in  calculating  mass  loading.  The 
representative  radius  is  used  with  the  relative  humidity  at  the  reference  level  to  solve  the 
modified  Hanel  droplet  growth  equation  iteratively  for  the  radius  of  the  dry  particles  of  this 
size  class.  Once  determined,  the  radius  of  dry  particles  for  each  size  class  remains 
constant  throughout  the  modeling.  However,  the  representative  droplet  radius  for  each 
size  class  of  evolving  droplets  is  allowed  to  change  without  regarding  ra  and  r^ 


4.9  Solution  Matrix 

The  system  of  equations  that  defines  the  model  is  written  and  solved  in  matrix  form.  At 
the  reference  height,  the  equations  consist  of  the  equilibrium  particle  growth  formulation 
(which  can  be  obtained  by  substituting  equations  (5),  (12),  and  (13)  into  equation  (7)  and 
differentiating);  equation  (44),  the  bulk  phase  change  relationship;  equation  (52),  the 
moisture  conservation  equation  including  the  drop  size  distribution;  and  equation  (22a), 
the  thermodynamics  of  the  ascending  sphere.  Above  the  reference  height,  the  equations 
consist  of  equation  (32),  the  nonequilibrium  growth  expression;  equation  (52);  and 
equation  (22a). 

To  obtain  reasonable  accuracy  in  reasonable  computation  time,  we  generally  divide  the 
drop  size  distribution  into  80  size  classes;  therefore,  the  basic  set  of  equations  consists 
of  82  equations  (83  at  the  reference  height).  For  high  accuracy,  1000  size  classes  are 
sometimes  used.  For  purposes  of  illustration,  a  matrix  equation  for  use  above  the  refer¬ 
ence  height  is  shown  below  for  a  compressed  set  of  equations  consisting  of  only  one  size 
class. 


1  A  B 

/  R  \  ^ 

dr,/dz 

C 

PwN1rr1  Pa  0 

0  r  **  1 

L  '•'p  J 

dWv/dz 

dT/dz 

= 

yn«+i 

YN,c+2 
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(55a) 


G,  H,  rtf,  and  rj  depend  on  the  droplet  size  and  thus  may  be  subscripted  to  indicate  the 
droplet  size  class.  At  the  reference  height,  all  droplets  are  assumed  to  be  in  equilibrium 
with  the  air.  At  each  level  above  the  reference  height,  the  growth  of  droplets  is  limited 
by  rates  of  diffusion  of  water  vapor  and  heat  transfer  to  and  from  droplet  surfaces. 

A,  B,  and  C  in  equations  (55a),  (55b),  and  (55c)  are  derived  in  a  way  that  assures  that 
droplet  growth  does  not  exceed  the  rate  allowed  by  diffusion.  G  and  H  in  equations  (55e) 
and  (55b)  are  values  (as  defined  just  after  equation  (21))  originating  from  the  modified 
Hanel  droplet  growth  model  (1981)  for  drops  that  are  in  equilibrium  with  the  atmosphere. 
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5.  MODEL  EVALUATION 


It  is  apparent  in  analyzing  our  model  equations  that  the  vertical  profiles  of  the  drop  size 
distribution  and  associated  parameters  can  vary  considerably,  but  we  feel  in  a  sensible 
way.  For  instance,  Shettle  and  Fenn  (1979)  require  that  their  drop  size  distributions, 
which  we  use  to  define  drop  size  distributions  for  our  model  at  the  reference  height,  vary 
significantly  as  a  function  of  relative  humidity,  visibility,  and  geography.  However,  we 
truncate  these  distributions  to  remove  the  smallest  drops  (those  that  would  have 
condensation  nuclei  with  equivalent  dry  sphere  radii  less  than  0.04  ^m)  because  of  scav¬ 
enging  effects  (Huser  and  Holloway,  1984;  HSnel,  1984).  Another  source  of  variability  is 
reported  by  Hanel  and  Lehmann  (1981),  who  show  that  the  equilibrium  sizes  of  drops  are 
sensitive  to  the  temperature  and  relative  humidity  of  their  immediate  environment  (which 
would  be  inside  the  sphere  in  the  present  study)  and  the  chemistry  of  the  dry  nuclei, 
which  also  vary  with  geography;  especially  significant  is  the  contrast  between  those 
continental  and  maritime  nuclei  that  belong  to  the  larger-radius  modes  of  bimodal  dis¬ 
tributions.  We,  in  turn,  have  shown  that  drops  in  a  nonequilibrium  state  with  their  imme¬ 
diate  environment  grow  quite  differently  than  do  those  assumed  to  be  in  an  equilibrium 
state.  Unfortunately,  some  results  in  the  literature  are  based  on  the  assumption  that 
drops  are  in  equilibrium  with  their  environment,  which  tends  to  exaggerate  the  size  of  the 
larger  droplets. 

Our  model  also  shows  that  the  growth  rates  of  the  drops  vary  according  to  the  ascent  rate 
of  the  sphere,  and  that  the  ascent  rate  of  the  sphere  varies  with  buoyancy,  mass  loading, 
entrainment,  and  the  initial  ascent  rate  of  the  sphere.  In  turn,  the  buoyancy  varies  with 
the  temperature  and  composition  of  the  environment;  the  mass  loading  varies  with  a 
number  of  factors,  including  stability  and  thermodynamics.  Finally,  the  bulk 
thermodynamics  is  a  function  of  temperature,  pressure,  and  relative  humidity.  Many  fac¬ 
tors  interact  in  a  very  complicated  way.  With  all  the  assumptions  and  interactions  that 
go  in  the  model,  it  is  gratifying  to  obtain  drop  size  distributions  that  make  any  sense  at 
all.  On  the  other  hand,  it  should  be  no  surprise  at  all  that  drop  size  distributions  vary  as 
much  as  reported  in  the  literature. 

When  we  started  to  evaluate  our  model,  we  were  surprised  to  find  that  no  (so  few) 
measurements  are  available  of  the  type  and  resolution  required.  The  largest  known  sets 
of  data  are  those  that  were  collected  in  Germany,  mentioned  in  the  Introduction,  and  they 
are  seriously  limited  for  several  reasons,  the  most  critical  being  that  reliable  temperature, 
relative  humidity,  and  wind  profiles  were  not  obtained  while  collecting  the  drop  size  and 
concentration  data.  Furthermore,  according  to  Lala  et  al.  (1987)  and  Loveland  and 
Lindberg  (1988),  the  drop  data  are  highly  suspect  in  the  very  low  and  high  LWC  ranges. 
Nevertheless,  we  examined  several  cases  of  measured  data  from  Meppen  with  emphasis 
on  vertical  profiles  of  liquid  water  content  and  drop  size  distributions. 

We  viewed  the  very  low  stratus  clouds  as  if  they  have  characteristics  lying  between  fog 
and  higher  stratiform  clouds  and,  in  some  cases,  even  low  cumulus,  for  which  there  are 
a  number  of  reported  cases  in  the  literature.  Our  comparisons  are  both  qualitative  and 
quantitative  where  possible.  We  chose  the  following  parameters  for  comparison:  tem¬ 
perature,  relative  humidity,  liquid  water  content,  and  drop  size  distribution. 

The  comparisons  that  follow  are  based  on  sets  of  realistic  conditions  that  we  specified  to 
initialize  our  model.  For  all  three  modeled  examples  used  here,  Shettle-Fenn  distribution 
parameters  for  a  maritime  air  mass  were  used,  the  initial  value  for  the  ascent  rate  of  the 
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cluster  sphere  was  3  cm/s,  and  the  sphere  radius  at  the  reference  (initial)  height  was  2 
m  (which  is  the  same  value  as  the  reference  height  above  ground  level).  Other  values 
used  as  input  to  the  model  are  listed  in  table  6. 

TABLE  6.  INPUT  FOR  RK  MICROPHYSICS  MODEL 


Corresponding  Meppen  80  profile  no.  --> 

Lower  cutoff  (fim)  * 

Entrainment  parameter  ae 
a  in  cluster  lapse  rate  equation 
b  in  cluster  lapse  rate  equation 
c  in  cluster  lapse  rate  equation 
Initial  height  (m) 

Height  of  top  of  cloud  (m) 

Input  values  at  initial  height: 

Visibility  (km) 

Relative  humidity  in  cluster 
Pressure  (millibars)  ** 

Temperature  in  cluster  (K) 

Temperature  of  ambient  environment  (K) 


5 

6 

7 

0.04 

0.04 

0.04 

0.03 

0.03 

0.03 

1.4  x  10"9 

2  x  10_# 

1.4  x  10-’ 

8  x  10* 

6  x  103 

1  x  10" 

2  x  10"7 

8  x  10'* 

1.5  x  10"' 

2.00 

2.00 

2.00 

280 

320 

290 

1.9 

1.4 

1.02 

0.93 

0.944 

0.96 

1016 

1015 

1015 

272.33 

272.36 

272.54 

272.25 

272.26 

272.24 

*  Lower  cutoff  is  the  smallest  dry  particle  (condensation  nucleus)  radius  ( nm )  allowed. 
**  Pressure  at  initial  height  is  the  total  pressure  (in  millibars)  due  to  air  and  water  vapor. 


5.1  Temperature  and  Relative  Humidity 

A  common  approach  for  computing  the  vertical  temperature  profile  below  a  cloud  (espe¬ 
cially  cumulus)  is  to  assume  that  from  the  surface  to  the  base  of  the  cloud  conditions  are 
dry  adiabatic,  and  in  the  cloud  conditions  are  moist  adiabatic.  Our  model  shows  that  the 
cluster  lapse  rate  is  close  to  dry  adiabatic  or  superadiabatic  at  reference  level  and  de¬ 
creases  to  approximately  moist  adiabatic  near  the  base  of  the  cloud.  In  the  lower  part 
of  the  cloud  where  relative  humidity  is  maximum,  the  lapse  rate  of  the  cluster  is  a  little 
greater  than  the  moist  adiabatic  value.  This  effect  and  the  superadiabatic  lapse  rate  at 
reference  level  both  appear  to  be  consequences  of  cooler,  drier  air  being  entrained  into 
the  cluster.  At  the  same  time,  our  model  shows  that  the  ambient  environment  lapse  rate 
may  vary  from  slightly  greater  than  dry  adiabatic  to  between  dry  adiabatic  and  moist 
adiabatic  in  the  subcloud  region.  The  cluster  and  ambient  atmosphere  (environmental) 
lapse  rates  are  shown  in  figure  5,  and  the  temperature  profiles  in  figure  6  for  the  three 
cases  that  we  considered  in  this  study.  Note  from  figure  6  that  it  would  be  difficult  to 
measure  the  differences  between  cluster  and  ambient  temperatures,  which  are  never 
greater  than  0.3  K  for  these  examples.  Noonkester  (1982)  presents  measured  temper¬ 
ature  and  relative  humidity  profiles  for  several  maritime  stratus  clouds.  We  chose  one 
case  from  his  figure  2  for  illustrative  purposes,  figure  7.  Note  that  the  vertical  structure 
of  both  the  temperature  and  relative  humidity  curves  is  generally  in  line  with  our  model 
results.  Also,  note  that  Noonkester's  data  were  collected  during  warmer  conditions,  so 
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that  the  moist  adiabatic  lapse  rate  is  smaller  than  some  of  our  cases  when  the  temper¬ 
ature  was  near  zero  degrees  centigrade.  For  stratocumulus  clouds,  Kao  and  Yamada 
(1989),  Nicholls  (1984),  and  Slingo  et  al.  (1982b)  found  that  the  temperature  gradient  is 
approximately  dry  adiabatic  in  the  subcloud  region  and  near  moist  adiabatic  in  the  cloud. 
For  stratocumulus  and  stratus  clouds,  Borovikov  (1963)  concluded  that  under  the  clouds 
the  lapse  rate  remained  lower  than  the  dry  adiabatic  lapse  rate  all  year.  Whereas,  in  the 
cloud  there  are  two  kinds:  either  a  more  or  less  uniform  lapse  rate  approaching  the  moist 
adiabatic  lapse  rate,  or  a  thin  inversion  layer  is  present  in  the  upper  region  of  the  cloud 
layer  and  the  cloud  "wedges”  into  the  overlying  main  inversion.  From  the  above  results 
we  concluded  that  the  temperature  profiles  provided  by  our  model  are  well  in  line  with 
what  should  be  expected. 

In  addition  to  the  temperature  profiles  we  also  computed  the  potential  temperature  gra¬ 
dient  profiles  of  the  cluster  and  ambient  atmosphere  for  the  cases  presented  in  figures  5 
and  6.  Results  are  shown  in  figure  8.  Note  that  the  potential  temperature  gradient  of  the 
cluster  is  negative  from  the  surface  to  approximately  77,  73,  and  167  m  above  ground  level 
for  profiles  5,  6,  and  7,  respectively,  and  then  becomes  positive,  suggesting  unstable 
conditions  in  the  lower  part  of  the  profile  and  stable  conditions  above  these  heights.  In 
the  ambient  environment,  the  corresponding  gradients  are  negative  from  the  surface  up 
to  about  34,  32,  and  64  m  above  ground  level  and  positive  above  these  heights,  suggest¬ 
ing  stable  conditions  start  at  lower  levels  in  the  ambient  environment. 

5.2  Liquid  Water  Content 

We  use  the  term  mass  loading  to  include  both  liquid  water  and  dry  particle  (nuclei)  ma¬ 
terial.  Field  measurements  are  actually  mass  loading  values.  In  our  model  we  can  sep¬ 
arate  the  contributions  of  each;  however,  the  field  measurements  do  not  lend  themselves 
to  this  treatment.  There  are  usually  no  serious  violations  if  one  assumes  that  the  drops 
are  all  water  in  the  truncated  distribution  used  in  our  model.  However,  caution  must  be 
used,  for  example,  when  modeling  absorption  for  wavelengths  at  which  the  imaginary  part 
of  the  index  of  refraction  is  much  larger  for  the  nuclear  material  than  for  pure  water. 

Estimates  of  liquid  water  content  in  clouds  are  most  controversial.  For  us  they  are  most 
frustrating.  One  can  defend  within  reason  any  estimate  of  LWC  in  clouds  that  one 
chooses  based  on  results  published  in  the  literature.  The  situation  is  even  worse  in  the 
subcloud  region.  Warner  (1970),  who  is  quoted  quite  frequently  in  the  literature  (for  in¬ 
stance  Pruppacher  and  Klett,  1978;  Rogers  et  at.,  1985);  references  studies  by  Ackerman 
(1959,  1963),  Skatski  (1965),  and  Warner  and  Squires  (1958)  from  which  he  states  that  the 
ratio  of  measured  to  adiabatic  liquid  water  decreases  with  height  above  cloud  base.  In 
contrast,  Driedonks  and  Ouynkerke  (1989),  referencing  Slingo  et  al.  (1982  a,b)  and  Nicholls 
(1984)  concludes  that  for  stratocumulus  clouds  the  mean  liquid  water  content  distribution 
is  very  close  to  the  adiabatic  value,  and  according  to  Nicholls  (1984)  earlier  measure¬ 
ments  of  liquid  water  content  that  gave  values  much  smaller  than  the  adiabatic  ones  are 
not  reliable  due  to  instrumental  inaccuracies.  However,  according  to  Slingo  et  al.  (1982b), 
who  studied  nocturnal  stratocumulus,  the  liquid  water  path  from  a  microwave  radiometer 
on  the  average  gave  values  0.7  times  the  value  for  adiabatic  ascent.  Rogers  et  al.  (1985) 
state  that  for  stratiform  clouds  (St,  Ns  and  As)  the  adiabatically  generated  LWC  is  used 
at  its  full  value  since  entrainment  is  relatively  unimportant  in  wide  spread  layered  clouds. 
At  the  other  end  of  the  spectrum  are  conclusions  by  Mao  Jie-tai  et  al.  (date  unknown), 
who  reference  Borovikov  (1961),  Singleton  and  Smith  (1960)  and  measurements  at  Wang- 
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jih-tai  during  the  period  September  1960  to  March  1961,  that  the  liquid  water  content  in 
layer  clouds  exceeded  adiabatic  values  by  as  much  as  a  factor  of  2. 

When  studying  cloud  condensation  nuclei  (CCN)  and  supersaturation  of  stratus  clouds, 
Hudson  (1983)  concluded  that  stratus  clouds  can  have  strong  adiabatic  characteristics. 

Liquid  water  content  estimates  from  recent  measurements  using  Knolfenberg  sensors  are 
supported  by  some,  adjusted  by  others,  and  argued  to  be  in  error  by  factors  of  2  to  4  in 
the  range  of  values  in  clouds  and  a  factor  of  5  in  subcloud  regions  by  Lala  (1987)  and 
Loveland  and  Lindberg  (1988). 

Low  (1980)  presents  quite  a  different  picture  of  the  LWC  profile  in  stratus  clouds  based 
on  data  from  the  Soviet  Union  (Khrgian,  1963)  and  data  from  Germany.  For  stratus  clouds 
defined  by  Koening  and  Schutz  (1974)  as  having  bases  ranging  from  0.1  to  1.0  km,  Low 
found  that  the  LWC  varies  parabolically  with  height  above  the  base  of  the  cloud  as 

LWC  =  A(0)  +  A(1)z  +  A(2)z* 

where 

A(0)  =  0.1193 
A(1)  -  1.99  x  10-3 
A(2)  =  -1.475  x  10-* 

LWC  is  in  g  m-3  and  z  in  meters. 

The  value  of  A(0)  would  be  suspect  by  many  researchers  who  are  inclined  to  assume  that 
the  LWC  at  cloud  base  is  zero. 

In  the  examples  calculated  by  our  model  for  this  report,  the  LWC  values  were 
subadiabatic  in  the  subcloud  region  having  a  profile  curvature  consistent  with  cloud 
measurements  made  in  Germany  (see  figure  1).  In  the  cloud  itself  the  model  produced 
values  ranging  from  0.5  to  0.75  of  adiabatic  to  agree  with  measured  values,  except  near 
the  cloud  base  where  relative  humidity  is  a  maximum.  See  figure  9.  LWC  profiles  re¬ 
ported  by  Noonkester  (1984)  for  maritime  stratus  are  similar  in  structure  to  those  simu¬ 
lated  by  our  model  (figure  10).  Our  strategy  was  to  adjust  the  cloud  LWC  values  so  that 
they  would  be  in  line  with  experimental  results  by  introducing  a  radiation  effect  (see 
section  4.3,  especially  equations  (22a),  (22d),  and  (22h)).  At  this  time  it  is  not  clear  as  to 
how  much  adjustment  should  be  made  in  general,  but  a  reasonable  mechanism  for  doing 
so  is  in  the  model.  It  also  allows  for  entrainment  effects  due  to  dry  air  insertion  from 
above  the  cloud. 

5.3  Supersaturation 

Hudson  (1983)  concluded  that  the  supersaturation  of  stratus  clouds  is  consistent  with  the 
concept  of  a  moderate  updraft  intermediate  between  values  determined  for  fog  and 
cumulus  clouds.  For  stratus  clouds  with  bases  between  300  and  600  m,  he  found  that  the 
effective  supersaturation  is  on  the  order  of  0.17  percent  which  is  less  than  the  more 
strongly  convective  cumulus  clouds  [for  example,  up  to  1  percent  (Twomey,  1959)]  but 
greater  than  that  of  the  less  active  fogs  [  ~0.10  percent  (Fitzgerald,  1978;  Hudson,  1980; 
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Meyer  et  al.  1980)].  More  recently  Hudson  (1986)  discusses  the  results  of  (our  field 
measurements.  Two  of  these  involved  stratus  clouds  “impacted”  on  California  coastal 
mountains,  while  the  other  two  sites  (1425  and  1520  m  elevation)  involved  fogs  in  Reno, 
Nevada.  The  California  sites  were  Henninger  Flats  (765  m  elevation)  near  Pasadena,  and 
San  Marcos  pass  (660  m  elevation)  near  Santa  Barbara.  Hudson  found  that  the  effective 
supersaturation  ranged  from  0.02  to  0.05  percent  at  Henninger  Flats,  and  >  0.1  percent 
at  San  Marcos  Pass,  but  the  Henninger  Flats  values  were  questionable  due  to  measure¬ 
ment  resolution  problems.  The  maximum  supersaturation  produced  by  our  model  for  the 
three  cases  considered  in  this  study  (see  figure  11)  was  0.043  percent,  which  is  smaller 
than  most  of  those  suggested  by  Hudson.  This  value  was  obtained  for  simulated  profile 
7;  the  other  two  profiles  did  not  reach  saturation.  These  small  maximum  relative  humidity 
values  are  attributed  to  the  very  hygroscopic  nature  of  the  condensation  nucleus  type 
(that  of  Hanel-Lehmann  sample  number  4)  chosen  for  all  three  model  examples. 

5.4  Drop  Size  Distribution  Comparison 

We  introduce  this  section  with  some  comments  about  the  basic  features  of  drop  size 
distributions  as  we  perceive  and  use  them.  First,  we  choose  a  reference  height  drop  size 
distribution  as  defined  by  Shettle  and  Fenn  (1979)  for  specified  values  of  relative  humidity 
and  visibility.  We  then  truncate  this  distribution  to  remove  the  smallest  drops  of  the  dis¬ 
tribution. 

We  presume  that  the  smallest  drops  are  removed  from  the  distribution  due  to  scavenging 
as  discussed  by  Husar  and  Holloway  (1984)  and  Hanel  (1984).  The  dry  particle  size  dis¬ 
tribution  was  determined  by  using  the  truncated  Shettle-Fenn  distribution  and  assuming 
equilibrium  at  reference  height.  The  damp  reference  height  drop  size  distribution  is 
shown  in  figure  12;  the  dry  particle  size  distribution  at  the  reference  height  is  shown  in 
figure  13.  In  effect,  for  the  cases  presented  in  this  report,  our  droplets  have  a  minimum 
dry  radius  of  0.04  n m  and  maximum  dry  radii  of  about  53,  49,  and  45  for  modeled 
profiles  5,  6,  and  7,  respectively. 

The  drops  in  this  truncated  distribution  grow  in  accordance  with  our  formulations,  given 
in  section  4.2,  as  they  are  lifted  within  the  rising  cluster  (see  figure  14).  The  results  are 
generally  consistent  with  Hudson  (1983)  who  states  that  these  subcloud  concentrations 
consist  mostly  of  enlarged  haze  particles  that  have  not  reached  their  critical  size,  but  the 
relative  humidity  is  high  enough  to  produce  detectable  haze  droplets.  As  the  drops  con¬ 
tinue  to  rise  above  the  height  of  maximum  relative  humidity,  some  drops  may  be  small 
enough  to  start  evaporating  after  an  initial  period  of  growth,  larger  drops  may  become 
activated,’  and  the  largest  drops  are  growing  although  still  subsaturated.  In  a  number  of 
model  results,  drops  in  some  of  the  smaller  size  classes  appear  to  first  be  activated  but 
then  shrink  to  become  haze  droplets.  This  behavior  looks  strange  but  is  consistent  with 
the  model.  The  model  does  not  allow  water  vapor  to  diffuse  from  a  region  of  lower  rela¬ 
tive  humidity  to  a  region  of  higher  relative  humidity,  even  when  the  latter  is  in  equilibrium 
with  the  surface  of  a  drop  having  a  radius  greater  than  the  radius  at  critical  supersatu¬ 
ration.  The  relative  humidity  of  air  in  a  cluster  in  a  lower  region  of  the  cloud  exceeds  the 
critical  supersaturation  value  for  drops  in  these  size  classes.  Thus,  although  showing 
some  lag  in  growth  (compared  with  equilibrium  drop  size)  from  diffusion  limitation  on 


’  A  condensation  nucleus  is  said  to  be  "activated"  when  the  droplet  formed  around  it  reaches  its  critical  sire 
(Rogers  and  Yau,  1 989) 
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drop  growth,  drops  in  these  size  classes  become  activated.  However,  closer  to  the  top 
of  the  cloud,  the  relative  humidity  of  air  in  the  cluster  decreases  with  height  and  becomes 
less  than  the  relative  humidity  at  the  surface  of  the  drop.  Thus,  some  water  evaporates 
from  the  drop,  even  if  its  radius  is  greater  than  the  critical  supersaturation  value,  and  the 
drop  becomes  deactivated. 

The  fast  growth  rate  of  many  of  the  drops  is  consistent  with  Korolev  et  al.  (unknown)  who 
found  that  the  condensation  nuclei  begin  to  take  on  significant  water  at  approximately  10 
to  15  m  below  the  cloud  base.  However,  our  results  are  in  conflict  with  Hudson's  (1983) 
proposition  that  the  largest  drops  are  the  first  to  become  activated.  As  the  droplets  con¬ 
tinue  to  rise  we  notice  that  the  smaller  droplets  undergo  some  evaporation  while  the 
larger  drops  continue  to  grow.  Interestingly  enough,  the  largest  drops  do  not  achieve 
their  critical  supersaturation  values  up  to  the  cloud  top,  but,  even  so,  continue  to  be  larger 
than  those  that  do.  Hence  additional  spectrum  broadening  is  noted  and  the  average  drop 
size  increases  (figures  14  and  15).  The  model  results  are  supported  by  several  studies, 
for  instance,  Singleton  et  al.  (1960)  in  studying  low  layer  clouds  found  that  there  is  a  dis¬ 
tinct  trend  for  spectra  to  become  broader  as  cloud  thickness  increases.  In  his  paper  of 
CCN  measurements  within  clouds  Hudson  (1984)  notes  that  Baker  et  al.  (1980)  proposed 
that  inhomogeneous  mixing  resulted  in  the  evolution  of  broadened  droplet  spectra. 
Hudson,  however,  found  that  homogeneous  mixing  is  sufficient  to  explain  this  broadening, 
and  our  results  support  this  view.  Pruppacher  and  Klett  (1978)  also  state  that  the  drop 
spectrum  becomes  broader  as  the  cloud  thicl  -less  increases.  Slingo  et  al.  (1982b)  in 
observing  marine  stratocumulus  found  that  the  mean  drop  radius  increased  steadily  to¬ 
wards  cloud  top,  and  that  there  was  a  systematic  increase  of  mode  radius  from  about  5 
/jm  at  cloud  base  to  11  /im  at  cloud  top.  These  values  should  vary  depending  on  the 
height  of  the  cloud  base  and  the  height  of  the  cloud  top.  Korolev  et  al.  (unknown)  when 
observing  stratiform  clouds  found  that  at  the  condensation  level,  water  containing  nuclei, 
depending  upon  their  degree  of  ascent,  grow  rapidly,  and  the  average  droplet  size  in¬ 
creases.  Similarly,  Slingo  et  al.  (1982a)  observed  in  nocturnal  stratocumulus  that  in  two 
thick  clouds  the  drop  size  data  showed  an  almost  monotonic  increase  of  mean  radius 
from  cloud  base  to  top.  The  mode  radius  increased  from  about  3  /im  at  cloud  base  to 
about  8  /im  at  cloud  top.  Nicholls  (1984)  in  studying  a  stratocumulus  cloud  over  the  North 
Sea  found  an  increase  in  the  mean  volume  radius  (rv)  of  the  drops  as  shown  in  figure  16. 
In  this  case  the  cloud  base  was  approximately  400  m.  Slingo  et  al.  (1982b)  in  studying 
marine  stratocumulus  found  that  the  mean  radius  varied  as  shown  in  figure  17  using  data 
collected  with  an  ASSP.  Goodman's  (1976)  summary  of  California  coastal  fog  and  stratus 
shows  a  consistent  fog  microstruc.ture  with  mean  diameters  increasing  with  height,  drop 
size  distributions  broadening  with  height,  and  liquid  water  content  increasing  with  height. 

If  one  chooses  to  define  the  number  of  cloud  drops  as  those  that  are  activated,  this 
number  is  zero  until  the  cluster  gets  above  the  height  of  maximum  relative  humidity, 
which  in  our  case  is  a  relatively  short  distance  above  cloud  base.  Above  the  height  of 
maximum  relative  humidity,  the  number  of  activated  drops  exhibits  a  few  fluctuations  of 
considerable  magnitude,  largely  owing  to  the  behavior  of  droplets  in  the  smallest  few  size 
classes  retained  in  the  distribution.  If  one  elects  to  define  the  number  as  a  sum  of  the 
activated  and  the  largest  drops,  the  number  is  not  much  different  from  the  number  of  ac¬ 
tivated  drops  since  there  are  so  few  in  the  largest  drop  sizes.  Nole,  however,  that  this 
count  is  sensitive  to  the  behavior  of  the  smallest  size  classes.  In  computing  extinction  for 
different  wavelengths,  we  must  also  be  concerned  with  the  smaller  droplets  even  though 
they  are  not  activated.  Finally,  the  measured  drop  size  distribution  is  strongly  a  function 
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of  the  sensors  used,  especially  since  there  are  large  numbers  of  drops  in  the  smallest 
size  classes. 

These  features  are  illustrated  in  figure  18  where  we  plotted  the  concentration  of  drops  as 
a  function  of  height,  resulting  from  our  model,  for  lower  cutoff  values  of  0.25,  0.5,  and  1.0 
Atm  for  the  drop  radius  for  profile  7.  As  a  general  qualitative  comparison  we  show  meas¬ 
ured  concentration  profiles.  For  instance,  in  figure  19  we  have  a  drop  concentration  pro¬ 
file  of  a  marine  stratocumulus  cloud,  observed  by  Slingo  et  al.  (1982b)  using  an  ASSP. 
The  base  of  the  cloud  was  approximately  500  m  and  the  top  was  900  m.  Similarly,  we 
show  drop  concentrations  (number  density)  for  three  cases  (figures  20,  21,  and  22)  of 
nocturnal  stratocumulus  studied  by  Slingo  et  al.  (1982a).  The  drops  were  measured  with 
an  ASSP  probe  with  a  cutoff  of  1  ^m  radius.  Nicholls  (1984)  presents  a  drop  concentration 
profile,  figure  23,  for  stratocumulus  using  an  ASSP  (cutoff  1  /tm).  Figures  24  and  25  show 
a  comparison  of  continental  (San  Diego)  and  maritime  droplet  concentrations  measured 
by  Hudson  (1983)  with  an  ASSP  (cutoff  0.3  nm  radius).  Note  that  the  measured  data  had 
lower  drop  cutoff  radii  as  indicated,  because  of  instrumental  range;  whereas  the  simu¬ 
lated  data  (such  as  shown  in  figure  18)  had  lower  cutoffs  of  whatever  wet  radii  corre¬ 
sponded  to  dry  radii  of  0.04  /tm.  These  smallest  values  of  wet  radii  vary  with  relative 
humidity,  which  in  turn  varies  with  height. 

Knollenberg  aerosol  spectrometers  of  the  type  (FSSP-100C)  used  to  acquire  Meppen  80 
data  (see  Lindberg,  1982)  are  known  to  underestimate  drop  concentrations  when  the  liq¬ 
uid  water  content  is  low.  According  to  Lala  et  al.  (1987),  and  independently  by  Loveland 
and  Lindberg  (1988),  the  liquid  water  content  derived  from  measured  drop  size  and  num¬ 
ber  data  can  be  in  error  by  as  much  as  a  factor  of  10  when  the  liquid  water  content  is  low. 
Hence,  we  concluded  that,  for  the  subcloud  region  in  cases  presented  in  this  study,  the 
liquid  water  content,  drop  number,  and  mean  volume  radius  derived  from  the  measured 
data  are  in  error,  likely  being  worse  near  the  ground  and  better  near  the  cloud  base. 
Furthermore,  the  visibility  estimates  would  also  be  in  error.  For  instance,  for  profile  5  the 
measured  visibility  at  the  lowest  height  was  1.9  km;  whereas  the  visibility  computed  by 
applying  Mie  theory  (to  spheres  of  pure  water  using  measured  drop  size  and  number 
data)  and  Koschmieder's  relationship  was  13.2  km.  In  contrast,  based  on  the  studies  by 
Lala  et  al.  (1987)  and  Loveland  and  Lindberg  (1988),  we  concluded  that  the  drop  size  and 
number,  and  consequently  liquid  water  content  and  mean  volume  radius,  estimates  are 
reasonable  in  the  cloud 


SUMMARY  AND  CONCLUSIONS 

We  feel  that  we  have  structured  a  theoretically  based  model  that  provides  reasonable 
internally  consistent  “average"  vertical  profiles  of  temperature,  pressure,  density,  poten¬ 
tial  temperature,  specific  humidity,  relative  humidity,  liquid  water  content,  and  drop  size 
distribution  in  very  low  stratus  clouds  and  their  associated  subcloud  regions.  We  also 
feel  that  combinations  of  these  profiles  can  be  useful  to  Army  users,  including  research¬ 
ers  concerned  with  moisture  effects  on  EO  sensors,  systems  and  electromagnetic  propa¬ 
gation.  Furthermore,  we  feel  that  we  achieved  the  goal  of  designing  a  model  that  requires 
a  minimum  number  of  conventional  surface  (2  m  above  ground)  meteorological  inputs, 
that  is,  temperature,  pressure,  and  relative  humidity.  However,  we  did  find  it  necessary 
to  know  (specify)  the  height  of  the  cloud  top;  type  of  nucleus;  and  the  rate  of  ascent,  ra¬ 
dius  of  cluster  sphere,  and  a  visibility  parameter  at  2  m  above  ground.  We  certainly  did 
not  achieve  the  simplification  of  input  attained  in  models  by  Duncan  (1980)  and  Heaps 
(1982)  who  require  only  one  surface  input  to  compute  "gross  average"  empirically  based 
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liquid  water  content  or  extinction  profiles.  However,  it  is  apparent  from  analyzing  our 
model  that  because  of  the  number  of  inputs  required  we  can  expect  significant  variation 
in  the  profiles,  particularly  on  a  worldwide  scope.  Lastly,  we  also  managed  to  structure 
the  model  so  that  it  can  be  executed  on  a  comparatively  small  computer  (HP9000-840). 

Most  of  the  vertical  profiles  produced  by  the  model  are  generally  in  line  with  field  obser¬ 
vations.  However,  for  some  cases  there  is  considerable  variability  in  the  observed  liquid 
water  content  profiles.  This  is  unfortunate,  since  we  had  decided  to  adjust  our  model  so 
that  the  model  liquid  water  content  profiles  would  generally  agree  with  field  observations. 
Furthermore,  there  is  such  a  large  variation  in  the  observations,  and  controversy  by  au¬ 
thors,  that  we  are  not  sure  of  the  best  way  to  adjust  the  liquid  water  content  profiles.  We 
have  used  selected  liquid  water  content  profiles  calculated  from  measured  data  to  sug¬ 
gest  desired  values  in  the  cloud  region. 

Our  vertical  profiles  of  drop  size  distribution  are  generally  in  line  with  the  literature.  We 
show  a  broadening  of  the  distribution  above  the  height  of  maximum  relative  humidity  and 
an  increase  in  the  average  or  mode  drop  radius  with  height.  It  is  important  to  note  that 
the  size  distribution  of  drops  is  also  changing  with  height. 

In  the  subcloud  region  all  droplets,  since  they  are  in  subsaturated  air,  are  damp  haze  type 
that  grow  in  size  monotonically  from  the  surface  to  the  base  of  the  cloud.  Once  in  the 
cloud,  the  smallest  of  these  haze  droplets  tend  to  evaporate  somewhat  so  that  the  drops 
decrease  in  size  with  height.  In  the  version  of  the  model  presented  in  this  report,  these 
smallest  droplets  were  among  those  omitted  from  the  truncated  distribution. 

The  following  features  and  resuits  of  our  model  are  different  from  some  in  the  literature. 
Our  model  is  a  nonequilibrium  model.  That  is  to  say,  the  smallest  drops  of  the  distribution 
are  very  close  to  being  in  equilibrium  with  their  environment  (in  cluster),  but  all  the  larger 
ones  are  in  nonequilibrium  to  some  degree:  Rates  of  drop  growth  are  limited  by  the  rates 
of  diffusion  of  water  vapor  to  the  drop  surfaces.  This  is  especially  significant  since  the 
largest  drops  are  not  allowed  to  grow  as  large  as  they  would  if  they  were  permitted  to  be 
in  equilibrium.  This  leads  us  to  another  point  of  difference.  Hudson  (1983)  has  concluded 
that  the  largest  drops  of  the  distribution  are  activated  first.  We  find,  in  contrast,  that  the 
midsized  drops  are  activated  first.  In  fact,  the  largest  drops  do  not  become  activated  by 
the  time  they  reach  the  top  of  the  cloud  in  our  examples.  Even  so,  they  are  always  larger 
in  size  than  the  activated  drops. 

Finally,  we  found  no  mathematical  or  physical  reasons  suggesting  that  the  vertical  pro¬ 
files  of  all  parameters  e  not  continuous  for  clouds  whose  bases  are  less  than  300  m. 
Hence,  we  feel  that  th  inification  of  the  cloud  and  subcloud  region  is  justified  as  was 
done  in  the  model  development. 
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Figure  1. 


LWC  profiles  from  observed  data  (Meppen,  Germany,  1980). 
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Figure  2.  Extinction  profiles  from  observed  data  for 

wavelengths  of  4  /<m  (top)  and  10.6  /im  (bottom). 
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Figure  3.  Relationship  between  the  extinction  coefficient  (0.55  pm) 
at  heights  z  and  z  +  20  m  (after  Duncan  et  a!.,  1980). 


Figure  4.  Relationship  between  liquid  water  content  at  heights 
z  and  z  +  20  m  (after  Duncan  et  al.,  1980). 
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Figure  5a.  Cluster  and  ambient  air  lapse  rates  for  model  profile  5. 

(Solid  line:  cluster;  short-dashed  line:  ambient  air; 
long-dashed  line:  height  of  inflection  point  of  log  (LWC)  profile.) 


Figure  5b.  Cluster  and  ambient  air  lapse  rates  for  model  profile  6. 

(Solid  line:  cluster;  short-dashed  line:  ambient  air; 
long-dashed  line:  height  of  inflection  point  of  log  (LWC)  profile.) 
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Figure  5c.  Cluster  and  ambient  air  lapse  rates  for  model  profile  7. 

(Solid  fine:  cluster;  short-dashed  line:  ambient  air; 
long-dashed  line:  height  of  inflection  point  of  log  (LWC)  profile.) 


Figure  6a.  Cluster  and  ambient  temperature  profiles  for  model  profile  5. 
(Solid  line:  cluster;  short-dashed  line:  ambient; 
long-dashed  line:  height  of  Inflection  point  of  log  (LWC)  profile.) 
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Figure  6b.  Cluster  and  ambient  temperature  proxies  for  model  profile  6. 
(Solid  line:  cluster;  short-dashed  line:  ambient; 
long-dashed  line:  height  of  inflection  point  of  log  (LWC)  profile.) 


TEMPERATURE,  ketvins 


Figure  6c.  Cluster  and  ambient  temperature  profiles  for  model  profile  7. 
(Solid  line:  cluster;  short-dashed  line:  ambient; 
long-dashed  line:  height  of  inflection  point  of  log  (LWC)  profile.) 
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Figure  7.  Ambient  temperature  and  relative  humidity 
profiles  of  marine  stratus  (Noonkester  1982). 
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Figure  8a.  Potential  temperature  profiles  for  model  profile  5. 
(Solid  tine:  cluster;  dashed  line:  ambient  air.) 


Figure  8b.  Potential  temperature  profiles  for  model  profile  6. 
(Solid  line:  cluster;  dashed  line:  ambient  air.) 


52 


H0GHT, 


- x - 1 - 1 - 1  i - » - j - i - 1 - 1 - r 

-10  0.0  10  2.0  3.0  4.0  5.0  6.0  7.0  8.0  9.0  10.0  110 

dGny/dz  *10"' 


Figure  8c.  Potential  temperature  profiles  for  model  profile  7. 
(Solid  line:  cluster;  dashed  line:  ambient  air.) 


Figure  9.  Model  mass  loading  curves  for  model  profiles. 

(Solid  line:  profile  5;  long  dashes:  profile  6;  short  dashes:  profile  7.) 


53 


Figure  10.  LWC  profiles  of  Noonkester. 


SUPERSATURATION  ♦10"4 


Figure  11.  Model  relative  humidity  profiles. 

(Solid  line:  profile  5;  long  dashes:  profile  6;  short  dashes:  profile  7.) 
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Figure  12a.  Reference  level  drop  size  distribution  (damp)  for  model  profile  5. 


DROP  RADUS,  urn 

Figure  12b.  Reference  level  drop  size  distribution  (damp)  for  model  profile  6. 
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Figure  12c.  Reference  level  drop  size  distribution  (damp)  for  model  profile  7. 


Figure  13a.  Reference  level  dry  particle  size  distribution  for  model  profile  5. 
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Figure  13b.  Reference  level  dry  particle  size  distribution  for  model  profile  6. 
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Figure  13c.  Reference  level  dry  particle  size  distribution  for  model  profile  7. 
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DROP  RADIUS,  micrometers 


Figure  14c.  Mode)  drop  size  growth  for  profile  7. 
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DROP  RADIUS,  micrometers 

Figure  15a.  Variation  of  average  drop  size  with  height  for  model  profile  5. 

(Solid  line:  median  radius;  short-dashed  line:  mean  radius; 
long-dashed  line:  mean  volume  radius.) 


60 


Figure  15b.  Variation  of  average  drop  size  with  height  for  model  profile  6. 

(Solid  line:  median  radius;  short-dashed  line:  mean  radius; 
long-dashed  line:  mean  volume  radius.) 


Figure  15c.  Variation  of  average  drop  size  with  height  for  model  profile  7. 

(Solid  line:  median  radius;  short-dashed  tine:  mean  radius; 
long-dashed  line:  mean  volume  radius.) 
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Figure  18.  Model  drop  concentrations  when  lower  cutoffs  are 

(right  to  left)  0.25  pm,  0.5  pm,  and  1.0  pm  for  profile  7. 
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Figure  19.  Drop  concentration  profile  (Slingo  et  al.  1982b). 
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Number  density  (cmr3) 

Figure  20.  Drop  concentration  (Slingo  et  al.  1982a). 
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Number  density  (cm-3) 

Figure  21.  Drop  concentration  (Slingo  et  al.  1982a). 
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Number  density  (cm-3) 

Figure  22.  Drop  concentration  (Slingo  et  al.  1982a). 


Concentration  {cm- 3) 

Figure  23.  Drop  concentration  profile  (Nicholis  1984). 
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Figure  24.  Drop  concentration  profile  (Hudson  1983). 


Figure  25.  Drop  concentration  profile  (Hudson  1983). 
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APPENDIX 


HYDROSTATIC  APPROXIMATION 


Assumptions: 

a.  Hydrostatic  approximation  is  based  on  TOTAL  pressure  (contributions  from  both  dry 
air  and  water  vapor)  in  the  ENVIRONMENT  (ambient  air  outside  of  parcel). 

b.  TOTAL  pressure  inside  parcel  equals  TOTAL  pressure  in  the  environment  at  the  same 
height. 

Thus,  if  P  is  defined  to  be  total  pressure, 

P  =  Pd  +  e  =  Pd '  +  e'  (A  -  1) 

p  -  p«(1+-TL)  -  p«'(1+-TL)’  (A-2) 


where  t  -  0.621972 


-^7  =  -9r>\o,a,  =  -9(r>'a  +  W  >'a)  =  _g^a(l+w;) 
The  ideal  gas  law  is  used  in  substituting  for  pa ' 


-gPd'O+w;) 

RdT' 


(A -3) 


(A -4) 


Substituting  equation  (A-3)  into  equation  (A-4)  yields 
dP  =  -gP  /  uw;  \ 
dz  RdT'  V  1  +  w Jr.  ) 

Dividing  both  sides  by  P,  multiplying  both  sides  by  dz,  and  integrating  produces 


In  Pj  —  In  P|_1 


RdT'  V  1  +  Wv  '/r  / 


(A -5) 


(A -6) 


where  Pj  is  the  total  pressure  at  level  i,  P,  „i  is  the  total  pressure  at  level  i-1,  and  Az  is  the 
difference  in  heights  of  the  two  levels  (z,  —  z,  _i). 

T'  and  Wv  '  are  functions  of  z  but  are  approximated  by  the  constants.  T'  and  Wv  ',  respec¬ 
tively,  to  facilitate  integration.  The  smallness  of  the  fractional  changes  in  T  and  in  the 
fraction  involving  Wv '  makes  this  a  reasonable  approximation. 


Equation  (A-6)  leads  quickly  to 


Pi  =  P,_,  exp 


i  +  vy 
1  +  wv  '/r. 


(A -7) 
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Sometimes  it  is  desirable  to  have  an  expression  to  replace  dP^/dz.  From  equation  (A-2), 


Pd(l  + 


) 


=  P 


(A -8) 


Wy  \  dPd  pd  dwv 

£  /  dz  e  dz 


dP 

dz 


(A -9) 


Using  equation  (A-5)  to  substitute  for  dP/dz  and  then  substituting  for  P  according  to 
equation  (A-8)  and  dividing  both  sides  by  (1  4-  Wv/c)  produces 


£ ^ 

dz 


-gPd  /  i+w' 
r<jT'  l  i  +  w'v/. 


(A -10) 


If  desired,  the  right-hand  side  of  equation  (A-5)  can  be  expressed  as  .  where  Tv' 

is  the  virtual  temperature  in  the  environment.  Similarly,  Tv '  and  T„ '  can  be  used  in 
equations  (A-7)  and  (A-10). 
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